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ABSTRACT 

We present the results of a GBT survey of NH3(1,1), (2,2), (3,3) lines towards 631 Bolocam Galactic 
Plane Survey (BGPS) sources at a range of Galactic longitudes in the inner Galaxy. We have detected 
the NH3(1,1) line towards 72% of our targets (456), demonstrating that the high column density 
features identified in the BGPS and other continuum surveys accurately predict the presence of dense 
gas. We have determined kinematic distances and resolved the distance ambiguity for all BGPS sources 
detected in NH3. The BGPS sources trace the locations of the Scutum and Sagittarius spiral arms, 
with the number of sources peaking between i?Gai~ 4 — 5 kpc. We measure the physical properties 
of each source and find that depending on the distance, BGPS sources are primarily clumps, with 
some cores and clouds. We have examined the physical properties as a function of Galactocentric 
distance, and find a mean gas kinetic temperature of 15.6 K, and that the NH3 column density and 
abundance decrease by nearly an order of magnitude between i?Gai ~ 3 — 11 kpc. Comparing sources 
at similar distances demonstrates that the physical properties are indistinguishable, which suggests 
a similarity in clump structure across the Galactic disk. We have also compared the BGPS sources 
to criteria for efRcient star formation presented independently by Heiderman et al. and Lada et al., 
and massive star formation presented by Kauffmann et al. 48% of our sample should be forming stars 
(including massive stars) with high efficiency, and 87% contain subregions that should be efficiently 
forming stars. Indeed, we find that 67% of the sample exhibit signs of star formation activity based 
on an association with a mid-IR source. 

Subject headings: ISM: clouds — ISM: dust — ISM: molecules — star: formation 



1. INTRODUCTION 

Massive stars (M > 8M0) play a significant role in 
the evolution of their immediate environment, as well as 
evolution on a galaxy-wide scale via their many feedback 
processes including strong UV radiation that creates HII 
regions, outflows and winds which inject energy back into 
the surrounding medium, and supernova at the end of 
their lives (Kennicutt 2005). Given the importance of 
massive stars in galactic evolution, it is imperative that 
we understand how and where massive stars form within 
the Milky Way in order to intrepret observations of other 
galaxies. 

Large sample sizes are beneficial in this task, and many 
large samples have been chosen in various ways. For ex- 
ample. Plume, Jaffe & Evans (1992) selected 179 H2O 
masers from Cesaroni et al. (1988) and followed up with 
observations of dense gas tracers (CS J — 7—6 and 
CO J = 3 — 2). This H2O maser sample was followed 
up with lower J transitions of CS (Plume et al. 1997; 
Shirley et al. 2003) as well as submillimeter dust contin- 
uum (Mueller et al. 2002). Shirley et al. (2003) thor- 
oughly characterized the physical properties based on 
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CS J = 5 — 4 maps of a subsaniple of 64 sources from 
Plume et al. (1992). They found the following median 
properties: R = 0.32 pc, Afvir= 920 Mq, E = 0.60 g 
cm~^. Plume et al. (1997) found a mean volume density 
of (log(n/cm~^))= 5.9 based on multiple transitions of 
CS. Recently, Wu et al. (2010) have characterized prop- 
erties for a subsample of 50 massive, dense clumps from 
Plume et al. (1992) based on HCN and CS line emis- 
sion. They found a range for the mean of each prop- 
erty based on the different molecular line tracers used: 
(M^i,)= 1300 - 5300 Mq, 0.29 - 1.09 g cm'^, 

and {n)— 3.2 x 10** — 2.5 x 10^ cm~^. Other samples 
were defined based on the presence of ultra-compact HII 
(UCHII) regions and IRAS colors. Wood & Churchwell 
(1989b) complied a catalog of massive star-forming re- 
gions from the IRAS All Sky Survey based on typical FIR 
colors of the cool dust surrounding UCHII regions (Wood 
& ChurchweU 1989a). Beuther et al. (2002b) compared 
1.2 mm continuum observations and CS molecular line 
studies toward 69 UCHII regions based on the sample 
defined by Sridharan et al. (2002). Beuther et al. deter- 
mined physical properties similar to those of Shirley et 
al. (2003). 

These initial studies of massive star formation were 
important first steps toward characterizing large sam- 
ples, although the methods employed in identifying the 
sources resulted in highly biased samples. Requiring the 
presence of H2O masers or UCHII regions biases the sam- 
ples to a particular stage in star formation. Both phe- 
nomena require a significant source of energy, and there- 
fore will bias the samples to later evoluntionary states 
when a protostar has already formed. Therefore, the 
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properties deterniined by these initial studies may not 
be characteristic of massive star forming regions at aU 
evolutionary stages because of these biases. 

The many large-scale, blind Galactic plane surveys re- 
cently completed and currently underway provide the 
data sets needed to study massive star formation in a 
Galactic context. In particular, three large scale surveys 
have recently done so: The Boston-University Five Col- 
lege Radio Astronomy Observatory ^"^CO Galactic Ring 
Survey (GRS; Jackson et al. 2006), the Red MSX Source 
(RMS; Hoare et al. 2004; Urquhart et al. 2008) Survey, 
and the Bolocam Galactic Plane Survey (BGPS; Aguirre 
et al 2011). 

Rathborne et al. (2009) identified 829 clouds and 6124 
clumps within the ^"^CO emission in the GRS. They com- 
pared properties of clouds inside and outside of the 5 
kpc molecular ring, a Galactic structure identified as an 
overdensity in molecular gas located at a Galactocentric 
radius of 5 kpc that contains the majority of star forma- 
tion in the Milky Way (e.g. Burton et al. 1975; Scoville 
& Solomon 1975; Cohen & Thaddeus 1977; Robinson et 
al. 1984; Clemens et al. 1988; Kolpak et al. 2002). They 
found that clouds within the 5 kpc molecular ring have 
warmer temperatures, higher column densities, larger ar- 
eas, and contained more clumps than clouds located out- 
side of the molecular ring. Roman-Duval et al. (2009) re- 
solved the distance ambiguity to 750 of the GRS clouds 
and found that their positions are consistent with a four 
arm model of the Galaxy and trace the positions of the 
Scutum-Crux and Perseus spiral arms. Roman-Duval et 
al. (2010) derived physical properties of 580 of the GRS 
clouds and demonstrated that the gas surface density 
peaks within the 5 kpc molecular ring. They found cloud 
sizes ranging from 0.1 to 40 pc, masses ranging from 10 
to 10^ Mq, a mean volume density of 230 ± 21 cm^'^ 
with the highest density being ~ 10'^ cm~'^, and a mean 
surface density of 144 ± 3 Mq pc~^. 

The RMS team has identified a candidate list of 
^2000 massive young stellar objects (MYSOs) and ultra- 
compact HII (UCHII) regions based on MSX and IRAS 
colors (Lumsden et al. 2002). Urquhart et al. (2011) 
have recently compared the radial velocities of the 
MYSOs to the clouds identified in the GRS and adopted 
the kinematic distances determined by Roman-Duval et 
al. (2009) thereby obtaining a complete sample of 196 
RMS sources with luminosities above their completeness 
limit of ~ 10^ L0. The clouds containing a RMS source 
were found to be larger, more massive, and more turbu- 
lent than the GRS clouds without a RMS source. The 
subsamplc of GRS clouds containing RMS sources are 
also well correlated with the Galactic spiral arm struc- 
ture. 

The BGPS (Aguirre et al. 2011) has detected 1.1 mm 
thermal dust emission from the dense regions closely as- 
sociated with star formation. Dunham et al. (2010; here- 
after refered to as DIO) found mean H2 column densi- 
ties of 1.2 X 10^^ cm~^ for BGPS sources located within 
the Gemini OBI Molecular cloud, and Schlingman et 
al. (2011) found a mean H2 column density of 6.8 x 10^^ 
cm-2 for BGPS sources detected in HCO+(3-2). In 
contrast to previous large-scale studies of massive star- 
forming regions, the BGPS is identifying high column 
density regions regardless of star formation signposts 



(such as II2O maser emission). By identifying sources 
based on dense gas and dust, the bias toward a particu- 
lar evolutionary stage is avoided. 

In this paper, we present the results of a targeted sur- 
vey of NH3 toward 631 BGPS sources in the first Galactic 
quadrant. The goals of this paper are to determine the 
location of a large sample of BGPS sources within the 
Galaxy, characterize the physical properties of a large 
relatively unbiased sample of star-forming regions, study 
their physical properties as a function of environment 
using Galactocentric radius as a proxy, and to charac- 
terize their star formation potential and mid-IR source 
content. Section [5] briefiy describes the BGPS and source 
selection, and Section|3]describes the NH3 set up and ob- 
servations. We present basic results of the NII3 survey in 
SectionlH In Section[5]we determine kinematic distances 
and break the kinematic distance ambiguity. Section [6] 
includes a description of each physical quantity derived 
from the BGPS and NII3 surveys. Section [7] explores 
the trends in physical properties as a function of Galac- 
tocentric radius, and Section [8] includes a discussion of 
the star formation activity, comparison to other large- 
scale studies, and a discussion of whether BGPS sources 
are forming massive stars. Finally, Section [9] provides a 
summary of this work. 

2. THE BOLOCAM 1.1 MM GALACTIC PLANE SURVEY 

The BGPS^ has observed approximately 170 square de- 
grees of the northern Galactic plane in 1.1 mm continuum 
emission using Bolocam at the Caltech Submillimeter 
Observatory'* (CSO). The survey area consists of a con- 
tinuous region spanning —10° < i < 90° and |5| < 0.5° as 
well as select regions of known star formation in the outer 
Galaxy: the Perseus tangent, the W3/4/5 region, and 
the Gemini OBI molecular cloud. The survey methods 
and data reduction are thoroughly described in Aguirre 
et al. (2011), and the source extraction algorithm and 
catalog are described in Rosolowsky et al. (2010). The 
BGPS has identified 8,358 continuum emission sources 
over the entire coverage area (Rosolowsky et al. 2010), 
most of which were previously unknown. Here we discuss 
details relevant to this work and refer the reader to these 
papers for further details. 

The effective FWHM beam size of the BGPS is 33", 
slightly larger than the nominal Bolocam beam at 1.1 
mm due to combining multiple observations of each field. 
The FWHM beam size corresponds to a solid angle of 
2.9 X 10~* sr, which is equivalent to a topliat func- 
tion with a 40" diameter {n = 2.95 x 10"^ sr). The 
BGPS catalog algorithm, bolocat, presents aperture fiux 
densities within 40", 80", and 120" diameter apertures 
representing the fiux density within one (5,^(40")), two 
(5*1,(80")), and three (5*^(120")) times the beam, respec- 
tively. Bolocat additionally provides an integrated fiux 
density (6*1, (int)), which is simply the sum of all pixels 
assigned to a given BGPS source and more accurately 
represents the fiux densities of irregularly shaped sources. 
The aperture-based flux densities of point sources require 
correction factors (1.46 for 5*^(40"), 1.04 for 5^(80"), and 
1.01 for 5,^(120")) to account for power falling outside of 

See |http://irsa.ipac.caltech.edu/data/BOLOCAM_GPS/ | 
* The Caltech Submillimeter Observatory is supported by the 
NSF. 
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Observed 1.1 mm Properties 



ID 






RA 


Doc 


— 

peak RA 


peak Dec 






— 

PA 


^obi 


5,^(120") 


Sviint) 


iiuiiibcr 


oourcG 


J ZUUU ] 








\ ) 


K ) 


K ) 


K ) 


(mJy) 


(mJy) 


1307 


G007.501+00.001 


18 02 30.6 


-22 28 00.3 


18 02 30.0 


-22 28 07.0 


23.4 


18.2 


173 


35.41 


399 (110) 


356 (94) 


1309 


G007.509+00.403 


18 01 00.2 


-22 16 03.4 


18 01 00.4 


-22 15 46.3 


27.1 


15.3 


154 


28.73 


550 (140) 


588 (110) 


1310 


G007.564-00.042 


18 02 46.4 


-22 26 12.5 


18 02 47.8 


-22 26 05.6 


28.2 


21.6 


34 


48.07 


629 (170) 


533 (170) 


1311 


G007.600-00.142 


18 03 14.7 


-22 27 07.6 


18 03 15.0 


-22 27 10.1 


25.5 


15.0 


142 


25.78 


717 (160) 


681 (120) 


1313 


G007.622-00.000 


18 02 46.8 


-22 21 46.5 


18 02 45.8 


-22 21 49.4 


18.7 


8.9 


105 


< 16.5 


494 (170) 


159 (90) 


1314 


G007.632-00.110 


18 03 12.8 


-22 24 28.3 


18 03 11.9 


-22 24 33.1 


29.7 


15.8 


40 


33.26 


1340 (190) 


1450 (150) 


1316 


G007.636-00.194 


18 03 31.5 


-22 26 45.8 


18 03 31.4 


-22 26 49.4 


18.0 


14.9 


77 


18.37 


617 (140) 


595 (94) 


1317 


G007.650-00.166 


18 03 27.6 


-22 25 13.3 


18 03 26.8 


-22 25 15.9 


24.3 


17.6 


35 


34.95 


735 (150) 


562 (120) 


1320 


G007.762+00.076 


18 02 48.1 


-22 12 09.8 


18 02 46.5 


-22 12 16.0 


29.1 


19.6 


15 


44.81 


584 (180) 


565 (170) 


1322 


G007.898-00.012 


18 03 23.9 


-22 07 42.2 


18 03 23.7 


-22 07 45.9 


29.1 


14.7 


173 


25.12 


615 (170) 


552 (130) 



Note. — The full tabic is available in the online journal. 

the aperture due to the sidelobes of the beam. S'^(int) 
does not require a correction since it is not aperture- 
based. 

In addition to the aperture corrections, a flux cahbra- 
tion factor of 1.5 ± 0.15 is required. This factor is based 
on a comparison of BGPS data with 1.2 mm data (Rath- 
borne et al. 2006; Motte et al. 2003, 2007; Matthews et 
al. 2009) acquired at different facihties, which showed 
that the BGPS flux densities were systematically lower 
(Aguirre et al. 2011). We apply the flux calibration factor 
in order to bring the BGPS flux densities into agreement 
with the 1.2 mm surveys. The aperture corrections and 
flux calibration factor are included in the flux densities 
presented here. 

Bolocat includes coordinates of both a geometric cen- 
troid and the peak 1.1 mm emission. The major {umaj) 
and minor {amin) axes of each source are calculated 
based on emission weighted moments, and the effective 
radius (i?) of each source is calculated as the geometric 
mean of the deconvolved major and minor axes: 



where abeam = Qfwhm /Vs\n2 = 14", Qfwhm = 33", 
and 7y is a factor relating the axis dispersions to the true 
size of the source, r] = 2.4 is adopted as the median value 
derived from measuring the observed major and minor 
axes compared with the true radius for a variety of sim- 
ulated emission profiles spanning a range of density dis- 
tributions, sizes relative to the beam and signal-to-noise 
ratios (Rosolowsky et al. 2010). Including r/ provides an 
effective radius that corresponds to the full extent of the 
emission detected by the BGPS. Specifically, R describes 
the radius over which Si, (int) was calculated. Due to the 
cleaning methods employed (PCA cleaning; see Aguirre 
et al. 2011), the BGPS resolves out uniform structure on 
size scales larger than 5.9'. Thus, the BGPS is sensitive 
to structures with diameters between 33" and 5.9'. 

In order to sample source properties across a large 
range of Galactocentric radius (i?Gai)j we have obtained 
NH3 and H2O spectral line observations toward 631 
BGPS sources within four distinct Galactic longitude 
ranges: 7.5° <£< 10.0°, 19.5° <£< 22.5°, 31.3° <£< 
34.5°, and 52.5° < £ < 55.5°. Sources near i = 33° were 
targeted because of the peak in mm continuum sources 
along that line of sight, which corresponds to the end of 
the Galactic bar and tangent of the Scutum spiral arm. 



The £ ^ 54° field was chosen to sample the tangent of the 
Sagittarius Arm. The £ ~ 9° and ^ - 20° fields were cho- 
sen to sample sources with small iicai • Table [T] presents 
the 631 BGPS sources that lie within our Galactic longi- 
tude ranges: BGPS source catalog number (column 1), 
source name (column 2), centroid RA and DEC (columns 
3 and 4), RA and DEC of the peak emission (columns 
5 and 6), major and minor axes (columns 7 and 8), po- 
sition angle north of £ = 0° (column 9), radius (column 
10), S'^(120") and S'^(int) (columns 11 and 12). 

We have chosen to use S,^{120") as well as S'^(int) 
due to their complementary nature. 5^(120") is a well- 
defined quantity that is easily compared with other data 
sets but may not represent the true source flux due to 
source crowding and/or irregularly shaped, large sources. 
S',y(int) is difficult to reproduce and compare with other 
data sets since it is a product of the BGPS source extrac- 
tion algorithm. However, S'^(int) will more accurately 
represent the source flux density. 

When choosing which aperture size to use, we consid- 
ered the distribution of source sizes as well as the dis- 
tribution of angular separation between BGPS sources. 
The distribution of angular source sizes is shown in Fig- 
ure [1] where the white histogram represents the angu- 
lar source size of all 631 BGPS sources, and the gray 
histogram denotes sources unresolved compared to the 
beam. The dotted lines mark the 3 aperture sizes in- 
cluded in bolocat. The 120" aperture is larger than most 
sources and will therefore underestimate the source flux 
in the smallest number of BGPS sources, although the 
number of sources larger than the 120" aperture is not 
negligible. Figure [2] shows the distribution of the an- 
gular distance to the nearest BGPS source. This dis- 
tribution favors the smaller aperture sizes since they will 
suffer the least from source contamination. We must bal- 
ance the effects of missing flux and including additional 
flux when choosing an aperture size. As demonstrated 
in Figure [21 217 BGPS sources have at least one source 
closer than 120". We have chosen to use the 120" aper- 
ture even though the flux will be overestimated in 217 
of the 631 BGPS sources within our longitude ranges. 
5*^(120") will, on average, overestimate S',y(int) by 10%; 
the mean ratio of S'i,(120") to S'^(int) is 1.09 with a stan- 
dard deviation of 0.50. Therefore, we present 5^(120") 
as a comparison point for other studies, while our analy- 
sis is based on S'j,(int), which will not suffer from the two 
problems discussed. We stress that S,y{120") will over- 
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TABLE 2 
NH3 Observation Dates and 
Atmospheric Opacities 
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Fig. 1. — Distribution of angular source sizes {R) for the 631 
BGPS sources that fall within the selected ranges (see section [2]l. 
The white histogram denotes all source sizes, the gray histogram 
denotes sources that are unresolved compared to the BGPS effec- 
tive beam size (FWHM of 33", R = 16.5"), and the dotted lines 
mark the radii of the apertures presented in the BGPS catalog. 
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Fig. 2. — Distribution of angular distance to the nearest BGPS 
source. The dotted lines represent the diameters of the apertures 
presented in the BGPS catalog. 

estimate the fluxes to approximately one third of BGPS 
sources included in this paper. 

We additionahy consider S^{4:0") as a representative 
of the mm flux per beam since the solid angle of the 
40" aperture (2.95 x 10~^ sr) is equivalent to the solid 
angle subtended by the effective Gaussian beam of the 
BGPS (2.9 X 10"* sr). S^{40") is only employed when 
calculating the H2 column density per beam and NH3 
abundance (see Section 1^771) . 

3. NH3 SURVEY OF BGPS SOURCES 

We have conducted spectroscopic follow-up observa- 
tions of the 631 BGPS sources selected as described in 



the previous section. We observed all sources in the low- 
est inversion transitions of NH3 and the 22 GHz water 
maser transition using the Robert F. Byrd Green Bank 
Telescope^ (GBT). The dates of observations and 23 GHz 
zenith opacities are listed in Table [2] The observations 
consisted of two minute integrations on each source, with 
pointing and focus checks approximately every 1.5 hrs. 
We additionally checked the flux calibration by period- 
ically injecting a noise signal through observations of 
3C309.1, an NRAO flux calibrator with a known flux. 
We repeated this flux calibration observation during all 
sessions except for sessions lasting two hours or less (ses- 
sions 4 and 11). We did not observe any significant vari- 
ation throughout the measurements (< 10%), and use a 
mean calibration for all observations. 

The GBT spectrometer served as the backend and was 
configured to observe 4 intermediate frequencies (IFs). 
The majority of sources were observed with the follow- 
ing four lines: 22 GHz H2O, NH3(1,1), NH3(2,2), and 
NH3(3,3). Thirty of the spectra were observed with an 
identical setup by the Red MSX Source team (RMS), 
who were kind enough to share their raw data. The 
sources observed during 2008 were observed with the 
NH3(4,4) in place of NH3(3,3). The frequencies of each 
line are listed in Table[31 We utilized frequency switching 
with a 5 MHz throw, and the RMS team used a throw 
of 8 MHz. At these frequencies, the GBT FWHM is 31" 
and is well matched to the effective BGPS FWHM beam 
size of 33". 

We reduced the spectral line data using the GBTIDL 
reduction package. We extracted each scan, and folded 
it to account for the effects of frequency switching. In or- 
der to place the data on the Ta scale, we calibrated the 
noise diodes using the observations of sources with known 
fluxes. We then subtracted a first order baseline (second 
order for the H2O spectra), and converted the spectra 
to the scale by accounting for the atmospheric opac- 
ity. We then divided by the GBT main beam efficiency, 
77 ^ 0.81 at 23 GHz, to place the spectra on the T^b 
scale. Each polarization was calibrated separately then 

^ The GBT is operated by the National Radio Astronomy Ob- 
servatory, which is a facility of the National Science Foundation, 
operated under cooperative agreement by Associated Universities, 
Inc. 
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TABLE 3 

Spectral Line Observations 



Line 


u 

(GHz) 


mean RMS 
(mK) 


Stdev RMS 
(niK) 


NH3(1,1) 


23.6944961 


98 


28 


NH3(2,2) 


23.7226328 


93 


27 


NH3(3,3) 


23.8701270 


96 


26 


NH3(4,4) 


24.1394160 


141 


33 


H2O 


22.2350801 


87 


32 



averaged together, as were multiple observations when 
applicable. The mean rms and standard deviation of the 
rms for each possible IF are listed in Table [S] 

Figure [3] shows the observed spectra for BGPS source 
1322, including NH3(1,1), (2,2), (3,3) and the 22 GHz 
H2O maser transition. The green line overplotted in each 
panel displays the model fit described in Section [3] 

The parameter estimation follows the philosophy de- 
veloped by Rosolowsky et al. (2008) and adapted for use 
in DIG. The method models the emission as a beam- 
filling slab of ammonia with a variable column density 
(-^NHs), kinetic temperature (Tk), excitation tempera- 
ture (Tex), Gaussian line width {avisr), and LSR velocity 
(z;lsr)- Previous versions of the pipeline used the opti- 
cal depth of the (1,1) line as a free parameter and then 
calculated the column density. This model reverses the 
direction of the calculation, assuming the molecules are 
in thermodynamic equilibrium using an ortho-to-para ra- 
tio of 1:1, which is the high temperature formation limit 
(Takano et al. 2002). Hence, the ammonia molecules are 
partitioned among the energy levels as 



Zn = 



l+Y. 2(2 J + 1) e^p^^^ '^lBJiJ+l) + iC-B)J^]+AEiJ,K,i) I 

for J = 3,6,9, = 0,1, (2) 



J,K,i 



Zp 



(2J + l)exp{-^iI^^(^ 



-(C-B)J^]+l\E(J,K^ 



J,K,i 



for J = l,2,4,5,...;i = 0, 1, (3) 

where Zq and Zp correspond to ortho and para NH3, 
respectively. Here, AE{J, K,i) is the energy difference 
due to the splitting of the symmetric and anti-symmetric 
states. AE{J, K,l) corresponds to the antisymmetric 
state, which is AE/k ~ 1.1 K above the symmetric state 
{AE{J, K, 0) =0). The column density of the molecules 
in the N^m{J: K, i) state is thus 7Vnh3-^o(>^, K, i)/{2Zo) 
and Nj<iB.3Zp{J, K,i)/ {2Zp), where the factor of 2 
in the denominator accounts for our assumption of 
ortho:para=l:l. From the column densities in the indi- 
vidual states, we calculate the optical depths in the indi- 
vidual transitions. From here, the method follows that of 
previous work: we use the optical depth, hyperfine struc- 
ture, the velocity information and the excitation condi- 
tions to model the individual spectra. The free param- 
eters are optimized using the MPFIT least-squares min- 
imization routine including parameter bounds (Mark- 
wardt 2009). 

These small changes do produce differences in the de- 
rived parameters, reflecting optimizing sightly different 
functions. In addition, the optimization weights the lines 
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Fig. 3. — Spectra toward BGPS source 1322. Upper panel: 
NH3(1,1). Lower panel from left to right: 22 GHz H2O maser tran- 
sition, NH3(2,2), and NH3(3,3). The green lines show the model fit 
to the observed spectra. Similar figures for the 514 NH3 pointings 
detected in the (1,1) line (see section |4)| are included as an online 
only figure. 



slightly differently under the different schemes. The vari- 
ation largely reflects the approximate nature of the as- 
sumed model, namely that of a homogeneous slab with 
uniform properties. While more sophisticated models 
could produce better flts, they would require additional 
information such as mapping to constrain the distribu- 
tion of the individual lines. As an example, consider 
the flt to an optically thin object which has a temper- 
ature that decreases with radius. In this case, the the 
model would underestimates the (3,3) line emission by 
a noticeable margin. This is because the temperature 
flt is largely driven by the bright (1,1) and (2,2) lines, 
which can be excited throughout the object. The (3,3) 
line could show signiflcant excitation in the center of the 
object but this emission would be underestimated since 
the temperature parameter is largely determined by the 
strong emission from the other two lines. Furthermore, 
variations in the ortho-to-para ratios could affect the flts; 
however we would require substantially more information 
to resolve these variations. The parameter fits represent 
emission- weighted averages of the parameters under con- 
sideration. 

The changes to the parameter estimation have resulted 
in approximately 10% changes in the NHs-based param- 
eters. For example, the mean Tk of the sources in the 
Gemini OBI molecular cloud have decreased by roughly 
2 K from the values presented in DIO. Because of the 
changes, we have reanalyzed the Gem OBI sources as 
well, and present the results from the newest pipeline in 
the comparison in this work. 

4. RESULTS 

Of the 737 NH3 pointings, 722 overlap with BGPS 
sources in the vl.O catalog, while the remaining 15 point- 
ings were selected by eye from an earlier BGPS cata- 
log. The 15 unmatched pointings are excluded from our 
analysis because they do not have corresponding 1.1 mm 
data. Eight of the 15 unmatched sources do have a de- 
tection in the NH3(1,1) line, suggesting that some of the 
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Fig. 4. — Distributions of 5iy(120") (top panel) and 5,y(int) (bot- 
tom panel). The white histograms include all BGPS sources ob- 
served in NH3, the gray histograms include the full sample, and 
the striped histograms include the Tk subsample. 



low intensity BGPS emission is likely real. We have ad- 
ditionally excluded ten of the 737 NH3 pointings due to 
radio frequency interference (RFI) in the NH3(2,2) IF. 
Six of the ten pointings with RFI were detected in the 
NIl3(l,l) transition, but we are unable to characterize 
the gas properties without a realistic upper limit for the 
(2,2) line. 

To detect a line, we require W > SAW and Tpk> 4cr, 
where W is the integrated intensity, AW is the error 
in integrated intensity, and Tpk is the peak brightness 
temperature of the line. We have detected the NIl3(l,l) 
line in 514 of the 707 NH3 pointings that were matched 
with a BGPS source and were not affected by RFI, for a 
detection rate of 73%. Since the positions of the observa- 
tions obtained in 2008 were determined by eye from early 
BGPS data products there are multiple NH3 pointings 
per BGPS source in some cases. The 514 NH3 pointings 
detected in the (1,1) line correspond to 456 unique BGPS 
sources (72% of the 631 BGPS sources observed). 

Table 2] lists the observed NII3 properties for the 514 
pointings detected in the (1,1) line, including BGPS 
source number (column 1), RA and DEC of the NII3 
pointing (columns 2 and 3), radial velocity (Vi^SR; col- 
umn 4), velocity dispersion {crvisr, column 5), and peak 
main beam temperature (Tmb) and integrated intensity 
(W) for the (1,1), (2,2), (3,3), and (4,4) NH3 transitions 
(columns 6-13). Lines that were observed but not de- 
tected are listed as upper-limits based on the observed 
rms (cr), and are given by Tnib< 4tT and W < 5 AW, 
where AW = aSV^/N is the error in the integrated in- 
tensity, SV is the channel width in velocity, and N is the 
number of channels over which the rms was calculated. 
These 514 spectra are included in an online version of 
Figure [H 

In the following analysis, we consider only the NH3 
pointings that lie closest to the 1.1 mm peaks, and refer 
to the BGPS sources and these pointings as the "full sam- 
ple". We additionally define two subsamples from the 
full sample: the "Tk subsample" and the "near subsam- 
ple". The Tk subsample consists of 199 BGPS sources 
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Fig. 5. — Observed peak T^^, of the NH3(1,1) line versus the flux 
within a beam, £'^(40"). Circles mark sources with a detection in 
both the (1,1) and (2,2) lines while crosses mark sources with only 
an upper limit in the (2,2) line. 



that were detected in the (2,2) line. Since sources with- 
out a (2,2) detection can only provide an upper limit for 
the gas kinetic temperature, the Tk subsample is used 
throughout this section and section [6] to characterize the 
properties of sources excluding sources that only have 
upper limits. The near subsample consists of 199 BGPS 
sources that were placed at the near kinematic distance 
(see section [5]), and is used in section [7] to assess the 
trends in source properties as a function of positon within 
the Galaxy. 

The distribution of 1.1 mm flux densities is shown 
in Figure 21 We present flux densities within both a 
well-defined aperture (S'i.(120"); top panel) and the in- 
tegrated fiux density, which is dependent on the reduc- 
tion and source extraction algorithms (S',y(int); bottom 
panel). The white histograms denote the distribution of 
all observed BGPS sources, the gray histograms denote 
the full sample, and the striped histograms denote the 
Tk subsample. The high flux density end of the distri- 
bution is well represented in both the full sample and the 
subsample, while a significant fraction of the faint BGPS 
sources are lacking (1,1) and (2,2) detections. 

Figure [5] plots Tnib(l,l) versus the 1.1 mm flux per 
beam, S'i,(40"), for the full sample. The Tk subset is 
shown as black circles, while sources without a detection 
in the (2,2) line are shown crosses. As expected, the (2,2) 
line was not detected in pointings toward the weaker 1.1 
mm sources with faint (1,1) lines. 

The observed radial velocities for the full sample are 
shown in Figure [HI The black circles denote the Tk sub- 
sample and the crosses mark BGPS sources without a 
(2,2) detection. The dotted line denotes the maximum 
radial velocity expected from Galactic circular rotation 
as a function of i given by fisr.max = 0(^o)(l — sin{£)), 
where 0(i?o) is the orbital velocity of the sun about the 
Galactic center. The observed radial velocities follow the 
distribution of ^^CO{J = 1-0) seen in the BU-FCRAO 
Galactic Ring Survey (GRS; Jackson et al. 2006). There 
are 24 sources that show some evidence of two velocity 
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Fig. 6. — Observed vi^j. for the full sample. Black circles mark 
the Tk subsample, while crosses mark sources without a (2,2) de- 
tection. The dotted line denotes the maximum radial velocity 
expected from circular rotation given by J^jsr mao; = &{Ro){^ — 
s'm{£)). Sources with radial velocities in excess of tijsr mai: a^re 
placed at the tangent distance. 

components. These sources are noted in Column 8 of 
Table El 

We also measure the Gaussian width of the (1,1) line, 
avisr, and plot the distribution of the full sample (gray 
histogram) and the Tk subsample (striped histogram) in 
FigureH Note that ayisr = FWHM/ VWn2. The dashed 
line denotes the spectral resolution of the NH3 observa- 
tions of 0.08 km s~^. The full sample can be charac- 
terized by ((Ty;sr)=0-76±0.49 km s"-'^, and the Tk sub- 
sample by ((Tv'isr)=0.81±0.53 km s^^. We observe line 
widths up to (Tyisr ^ 4 km s~^. Of the 11 sources with 
o'visr ^ 2 km s~^, six are weak detections with poor 
model fits, one has two velocity components, and four 
truly do exhibit line widths greater than 2 km s^^. See 
Table|H]for statistical properties of (Tyisr for both the full 
sample and the Tk subsample. 

We defer the analysis of the H2 O maser observations to 
a later paper, but we do note the presence or absence of 
a water maser in Column 7 of Table [S] We have detected 
the 22 GHz H2O line towards 182 BGPS sources (40%). 

5. RESOLVING THE KINEMATIC DISTANCE AMBIGUITY 

Without distances to the BGPS sources, the observed 
properties can not be translated into physical proper- 
ties. The observed NH3 radial velocity in conjunction 
with the Galactic rotation model of Reid et al. (2009; 
Ro = 8.4 ± 0.6 kpc, 60 = 254 ± 16 km s'^) provides 
kinematic distances to all sources with a detection in 
the NH3(1,1) transition. However, the application of the 
Galactic rotation model introduces a distance ambiguity 
(known as the kinematic distance ambiguity, or KDA) 
for sources within the solar circle. A single radial veloc- 
ity provides two distances (refered to as a near and far 
distance) found on either side of the tangent point. Only 
sources found at the tangent velocity do not suffer from 
the KDA. We break the KDA based on the presence of in- 
frared dark clouds (IRDCs) or HI self- absorption (HISA). 

IRDCs are regions seen in absorption against the dif- 
fuse infrared background, some of which have been shown 
to have high densities (> 10^ cm"'^), cold temperatures 
(< 25 K), and high column densities (~ xq^'^"^''^ cm~^; 
Egan et al. 1998; Carey et al. 1998, 2000; Simon et 
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Fig. 7. — Distribution of Gaussian line widths for the full sample 
(gray histogram) and the Tk subsample (striped histogram) . Note 
(^Vlsr = FWHM/\/8in2. The dashed line denotes the resolution 
of the NH3 observations of 0.08 km s~^. Six of the eleven sources 
with crvisr ^ 2 km s~^ are weak detections, one has two velocity 
components that are poorly fit with the current model, and the 
remaining four sources truly do exhibit line widths greater than 2 
km s^-"^ . 



al. 2006; Rathborne et al. 2005, 2006, 2007, 2008, 2010). 
Thousands have been cataloged from the MSX (Price et 
al. 2001) and GLIMPSE (Benjamin et al. 2003; Cham- 
bers et al. 2009) surveys (Egan et al. 1998; Carey et 
al. 1998; Simon et al. 2006; Peretto & Fuller 2009; Cham- 
bers et al. 2009). Spatial coincidence and structural sim- 
ilarity between a BGPS source and an IRDC places the 
BGPS source at the near kinematic distance where it 
is able to absorb the diffuse IR background. A BGPS 
source not coincident with an IRDC is assumed to be 
at the far kinematic distance. We have cross matched 
the BGPS sources observed in NH3 with the catalog of 
IRDCs seen in the Spitzer GLIMPSE (Peretto & Fuller 
2009) and identified 105 BGPS sources coincident with 
an IRDC. We have visually inspected the GLIMPSE 8 
fim images for spatial coincidence and structural similar- 
ity between the 1.1 mm emission and 8 /im absorption 
features both identified by Peretto & Fuller (2009) as 
well as those not included in their IRDC catalog. 47 
IRDCs from Peretto & Fuller (2009) that were matched 
to a BGPS source were deemed not to be evidence of the 
near distance due to a mismatch in position or structure 
of the IRDC compared to the BGPS emission. Over- 
all, 171 of the 456 BGPS sources in the full sample were 
found to coincide with an 8 fim IRDC and were placed 
at the near kinematic distance, while 215 were assigned 
the far distance due to a lack of correspondence between 
8 /im absorption and 1.1 mm emission. Additionally, 
70 BGPS sources had inconclusive evidence of an 8 
IRDC and were flagged as possibly at the near distance. 

We note that the assumption that an IRDC is located 
at the near distance could be false in some cases. For 
example, since the identification of an IRDC relies on 
a bright infrared background, a source at the near dis- 
tance without a significant infrared background to absorb 
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would be falsely assigned to the far kinematic distance. 
Similarly, sources exhibiting weak 8/zm absorption could 
be located at the far distance. Battersby et al. (2011) 
have identified a sample of 5 high column density sources 
in the Herschel Hi-GAL data that have weak 8/im ab- 
sorption and exhibit a discrepancy between the far-IR 
column density and that calculated from the 8 /im ab- 
sorption under the assumption of the near kinematic dis- 
tance, which suggests the sources may be located at the 
far distance. While our assumption that a BGPS source 
associated with an IRDC is at the near distance is ac- 
curate in most cases, it disregards the fact that some 
IRDCs may be located at the far distance. 

To minimize false kinematic distance assignments, we 
have also considered HI self-absorption when resolving 
the KDA. The HI self-absorption technique for resolving 
the KDA is similar to the IRDC method in that cold, 
dense material located at the near distance will absorb 
the emission from warmer material with the same radial 
velocity found at the far kinematic distance (e.g., Knapp 
1974; Burton et al. 1978). The signature of HIS A is fre- 
quently displayed in molecular clouds (e.g., Knapp 1974), 
and theoretical studies by Flynn et al. (2004) demon- 
strate that all molecular clouds contain enough HI to pro- 
duce a self-absorption profile against a warm continuum 
background. Indeed, HISA has been seen throughout the 
Canadian Galactic Plane Survey of 21 cm emission to- 
ward the outer Galaxy (Gibson et al. 2000, 2005). This 
technique has recently been used to resolve the KDA to- 
ward star-forming regions and HII regions (Jackson et 
al. 2002; Fish et al. 2003; Busfield et al. 2006; Pandian 
et al. 2009; Anderson & Bania 2009; Roman-Duval et 
al. 2010). There are situations in which this method 
for KDA resolution could be incorrect. For example, a 
BGPS source at the far distance could have a peculiar 
velocity that places it at the same velocity as a warmer 
source at the far distance and causes it to exhibit HISA. 
Also, an HII region embedded within a cold source at the 
far distance could also be incorrectly assigned to the near 
distance (Roman-Duval et al. 2010). In such cases, the 
BGPS source would be incorrectly assigned to the near 
distance. Here we accept the KDA resolution as correct, 
but note that it may not be correct in all cases. 

We have used the 21 cm VLA Galactic Plane Survey 
(VGPS; Stil et al. 2006) and Southern Galactic Plane 
Survey (SGPS; McClure-Griffiths et al. 2005) to search 
for HI self- absorption toward the 456 BGPS sources in 
the full sample. Figure [5] shows examples of possible 
outcomes of the HISA method. We have marked the ra- 
dial velocity observed in our NH3 observations, but use 
i^CO (1-0) spectra from the BU-FCRAO Galactic Ring 
Survey (GRS; Jackson et al. 2006) for comparison since 
^■^CO emission traces the line width and spatial extent of 
HISA well (e.g.. Goldsmith & Li 2003; Li & Goldsmith 
2005). The top panel shows a clear example of a HISA 
for BGPS source 2630. The green dashed line is the HI 
spectrum, the blue solid line is the ^'^CO spectrum, and 
the dotted black line marks the observed NH3 velocity. 
The HISA feature matches both the observed NH3 ve- 
locity and the ^^CO peak and line width. The middle 
panel shows a spectrum toward a BGPS source with po- 
tential HISA. The HI spectrum has an absorption feature 
slightly offset from the observed NH3 velocity and ^^CO 
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Fig. 8. — HISA examples. The green dashed line shows the HI 
spectrum, the blue solid line shows the ^'^CO spectrum (top and 
middle panels) or the observed NH3(1,1) spectrum (bottom panel), 
and the dotted black line marks the observed NH3 velocity. Top 
panel: BGPS source 2630. There is a very clear HISA feature at 
the same velocity as the ^-^00 peak and the observed NH3 velocity. 
Middle panel: BGPS source 2843 with potential HISA. The width 
of the potential HISA feature and the offset in velocity from the 
dense gas tracers resulted in a "m" flag for HISA. Bottom panel: 
BGPS source 1405 with no HISA. '^^CO observations were not 
available in the 7.5° < I < 10.0° field, we overlay our NH3 spec- 
trum instead. The two NH3 features seen in absorption on either 
side of the NH3 emission lines are due to the frequency switching 
used during observations. There is a clear lack of HISA at the same 
velocity as the NH3 spectrum. 



peak and the ^^CO line width does not match the poten- 
tial absorption feature. The bottom panel shows a clear 
example a BGPS source placed at the far distance based 
on the HISA method. Here the blue solid line displays 
the observed NH3(1,1) spectrum since the GRS did not 
include our 7.5° < I < 10.0° field. The NH3 hyperfine 
structure can clearly be seen in this spectrum, and the 
absorption features seen on either side of the NH3 emis- 
sion are a result of the frequency switching performed 
during observations. The HI spectrum is smoothly de- 
creasing at the NH3 velocity, and shows no absorption 
feature. Of the 456 BGPS sources in the full sample, 
the HISA technique has identified 116 BGPS sources at 
the near distance, while 130 displayed questionable HISA 
and could not definitively break the KDA and 210 sug- 
gested the far distance due to a lack of HISA. 
When investigating matches between a BGPS source 
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Fig. 9. — Distribution of kinematic distances. The gray his- 
togram marks the full sample of 456 BGPS sources, while the 
striped histogram denotes the Tk subsample. The peak in the 
distribution around 4 kpc is a result of the 5 kpc molecular ring. 



and either an IRDC or HISA we assign one of three flags 
to each source for each method: "y" denoting a positive 
association suggesting the near distance, "n" denoting 
no association, and "m" representing a questionable as- 
sociation. A BGPS source is assigned the near distance 
if it receives a single "y" (143 sources), two "y" flags (72 
sources), or two "m" flags, one from the IRDC and HISA 
methods (19 sources). Conversely, a BGPS source is as- 
signed the far distance if it receives two "n" flags (126 
sources) or one "n" and one "m" flag (96 sources). 

By combining the IRDC and HISA results, we have 
identified 233 BGPS sources at the near distance, 221 at 
the far distance, and two at the tangent distance where 
the KDA does not exist. Table [6] lists the BGPS source 
ID (column 1), HISA, IRDC, and KDA flags (columns 
2-4, respectively), as well as the assumed kinematic dis- 
tance (column 5; asymmetrical errors listed in paren- 
theses) and resulting i?Gai (column 6). Figure [3] shows 
the distribution of distances for the full sample of 456 
BGPS sources (gray histogram) and for the Tk subset 
(striped histogram). The peak in the distribution near 

4 kpc is likely a result of abundant dust and gas in the 

5 kpc molecular ring (e.g. Burton et al. 1975; Scoville 

6 Solomon 1975; Cohen & Thaddeus 1977; Robinson et 
al. 1984; Clemens et al. 1988; Kolpak et al. 2002; Rath- 
borne et al. 2009). Figure [TOl shows the location of the 
456 BGPS sources in the Galaxy. The green crosses mark 
sources placed at the near distance, red squares mark 
sources at the tangent distance, while blue triangles mark 
sources placed at the far distance. The dotted line de- 
notes the location of the tangent distance as a function 
of Galactic longitude. 

As previously discussed, the use of IRDCs and HISA 
to resolve the KDA may lead to incorrect assignments 
of both near and far kinematic distances. To assess the 
possible effects of incorrect distance assignments, we have 
explored the worst case scenario, in which all sources 
have been assigned incorrect distances, by reversing the 
KDA resolution for all sources in the full sample. This 
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Fig. 10. — Location of the full sample of 456 BGPS sources in the 
Galaxy. Green crosses mark BGPS sources placed at the near kine- 
matic distance, red squares mark sources at the tangent distance, 
and blue triangles mark sources placed at the far kinematic dis- 
tance. The dotted line marks the location of the tangent distance 
as a function of Galactic longitude. 



test showed little effect on the statistical properties of the 
sample, suggesting that our results will not be greatly 
affected by incorrect distance assignments. 

The use of a circular rotation curve for the Galaxy 
necessarily simplifies Galactic kinematics, and thus in- 
troduces errors into our assumed distances. For exam- 
ple, spiral arm structure can cause the radial velocities of 
nearby sources to deviate from circular rotation and dis- 
play peculiar velocities. Similarly, the radial velocities of 
nearby sources that are orbiting the Galaxy at the angu- 
lar velocities similar to that of the sun are dominated by 
peculiar velocities. Additionally, sources located within 
the Galactic bar exhibit streaming motions that can de- 
viate significantly from circular rotation (e.g. Burton & 
Bania 1974). In these cases, the observed radial ve- 
locities will return inaccurate kinematic distances un- 
der the assumption of simple circular rotation. Thus, 
kinematic distances have large inherent errors. Sources 
located in the £ ~ 8.5° field could be strongly affected 
by the streaming motions within the bar. Sources in 
this region assigned small kinematic distances should be 
viewed skeptically. It is possible that peculiar velocities 
or streaming motions have resulted in a very low radial 
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Fig . 11. — Distribution of Galactocentric radii for the full sample 
(gray histogram) and the Tk subsample (striped histogram). The 
peak in number of sources at i?Gal~4.5 kpc corresponds to sources 
in the 5 kpc molecular ring. 

velocity and the source is truly located within the Galac- 
tic bar or Scutum spiral arm at approximately 4-5 kpc. 
These uncertainities should be noted; however, without 
a detailed model of the spiral arm structure, peculiar ve- 
locities and streaming motions throughout the Galaxy, 
we are unable to correct for these effects. In this paper, 
we accept the kinematic distances as correct but note the 
potential for large errors. 

While the kinematic distance is subject to the near/far 
distance ambiguity in the inner Galaxy, each near/far 
distance pair corresponds to a single distance from the 
Galactic Center, the Galactocentric radius. Figure [H] 
plots the distribution of Galactocentric radius for the full 
sample (gray histogram) and the Tk subsample (striped 
histogram). The peak in the distribution at i?Gai= 4 — 
5 kpc corresponds to the 5 kpc molecular ring, which 
corresponds to the end of the Galactic bar and beginning 
of the Scutum spiral arm. The peak in source number 
at i?Gai=6 kpc corresponds to the Sagittarius arm. We 
do not have enough sources at large i?Gai for a definitive 
detection or non-detection of the Perseus arm at i?Gai=8 
kpc. Section[7]discusses trends in properties as a function 
of Galactocentric radius. 

6. ANALYSIS 

With kinematic distances, we are able to calculate 
physical properties such as mass and volume-averaged 
density. In this section, we present each property and 
describe the biases inherent in each. Errors are propa- 
gated through each derived quantity and are presented 
for individual sources within parentheses in Tables (TJ and 

HHII 

When discussing mean properties here and throughout, 
the uncertainties given are the standard deviation about 
the mean rather than the uncertainty in the mean. Some 
property distributions are non-Gaussian, and the uncer- 
tainty given should be viewed as a measure of the scatter 
in the sample rather than the error. 

General statistics (minimum, mean, standard devia- 



tion, median, and maximum) for the properties presented 
in this section for the full sample and the Tk subsample 
can be found in Tabled 

6.1. Physical Size 

We calculate the physical radius based on the object 
radius determined by the source extraction algorithm 
(Equation [T|) and the kinematic distance. The physical 
radii for the full sample are given in column 2 of Table 
[7| and are shown in Figure [T2k . The gray circles are the 
Tk subsample and the gray crosses are BGPS sources 
with no (2,2) detection. The black squares and the solid 
black line mark the mean physical radius within 2 kpc 
bins. The short dashed line marks the physical radii cor- 
responding to the beam radius of 16.5" as a function of 
distance, and the long dashed line marks the physical 
radii corresponding to 2.95' (half of the spatial scale at 
which the BGPS loses sensitivity due to spatial filter- 
ing). BGPS sources with extracted radii smaller than 
half the beam size have been set as upper limits with 
radii equal to half the beam size. The mean physical ra- 
dius increases with distance as expected due to the upper 
limit on spatial scales imposed by spatial filtering during 
data reduction and the lower limit due to the beam size. 

Typical radii limits for cores, clumps, and clouds are 
denoted with the dotted lines at 0.125 pc and 1.25 pc 
(cores R = 0.01 — 0.1 pc, clumps R — 0.15—1.5 pc, clouds 
i? = 2.0 - 7.5 pc; Bergin & Tafalla 2007) in Figure [H 
Assuming that BGPS sources have a constant density, we 
can place contraints on the distances at which each type 
of source would be included in the BGPS catalog. Cores 
farther than 1.25 kpc will be unresolved by the 33" beam 
(for R = 0.1 pc), although unresolved cores can still be 
detected by the BGPS if they have flux densities high 
enough to counter balance beam-averaging effects and 
remain above the beam-averaged detection threshold (see 
Section 7.2 of DIO). The largest {R = 1.5 pc) uniform 
density clumps at distances less than 1.75 kpc would be 
filtered out but would be resolved across the Galaxy. The 
smallest clumps {R — 0.15 pc) would become unresolved 
relative to the beam at 1.9 kpc. Similarly, the smallest 
clouds would be filtered out at distances closer than 1.2 
kpc, while the largest would be filtered out if closer than 
8.7 kpc. The largest clouds are resolved even at the far 
side of the Galaxy. 

The mean 1.1 mm source radius increases with distance 
and sources are extended compared to the 33" beam. 
There are 24 unresolved BGPS sources for which we will 
also need to consider mass in order to place them in the 
hierarchical structure. 

6.2. Gas Kinetic Temperature 

The gas kinetic temperature is calculated based on the 
excitation temperatures of the observed NH3 lines. For 
sources with only the (1,1) fine detected, the rms of the 
(2,2) spectra was used as an upper limit to calculate the 
(2,2) excitation temperature and the resulting Tk is only 
an upper limit. The uncertainties in individual Tk mea- 
surements are small, with (cj'j^/Tk) = 0.06. 

The derived Tk are given in column 3 of Table [51 The 
full sample is characterized by (Tk) < 15.6±5.0 K, while 
the Tk subsample is characterized by (Tk) = 17.4 ± 5.5 
K. The fractional scatter about the mean (0.32) is five 
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Fig. 12. — Physical properties versus kinematic distance. In all panels, the dotted lines delineate typical properties of clouds, clumps, and 
cores while the short dashed lines mark lower limits of the physical properties within a 33" beam. The gray circles mark the Tk subsample 
and gray crosses mark sources with no NH3(2,2) detection. The black squares and lines denote the mean of each property within 2 kpc 
wide bins centered at 1, 3, 5, 7, 9, 11, and 13 kpc. (a) Physical radius. The long dashed line denotes the upper limit to source size imposed 
by the spatial filtering of uniform emission on size scales greater than 5.9'. (b) -Miso(int). (c) Volume averaged density, (d) H2 column 
density. The dotted line in panel (d) marks the surface density of 1 g cm~^ required to prevent fragmentation and allow a massive star to 
form (Krumholz & McKee 2008). The dash-dot line denotes the threshold gas surface density required for efficient star formation as seen 
in nearby low mass star-forming regions (Lada et al. 2010; Ifciderman et al. 2010). 



times larger than the individual uncertainties in the mea- 
surements (0.06) so the differences among the sources 
are real. See Table [5] for further statistics. There are 
81 sources that likely require multiple Tk components to 
fit the NH3(3,3) emission as well as the (1,1) and (2,2). 
These sources are marked in Column 8 of Table [5l 

Without mapping the NH3 emission towards the BGPS 
sources we can not characterize kinetic temperature gra- 
dients in the gas. However, since the NH3 pointings 
correspond to the peak 1.1 mm emission, which closely 
traces the NH3 emission (Friesen et al. 2009), the mea- 
sured Tk is likely an upper limit as Tk has been shown to 
increase toward the peak of the NH3 emission (Zinchenko 
et al. 1997). Thus the NH3 outside of the GBT beam are 
likely to be colder than the Tk observed toward the peak 
of the 1.1 mm emission. 

6.3. Isothermal and Virial Masses 

Assuming the dust can be characterized by a single 
temperature, we calculate the isothermal mass, Miso, as 



= 13.1 (^) (r^)^( /r.%?:V, )(4) 

where Si, is the 1.1 mm flux density, D is the kinematic 
distance, Ki, is the dust opacity per gram of dust and gas 
and includes a gas-to-dust ratio of 100, and B,y is the 



Planck function evaluated at Td. We logarithmically in- 
terpolate the Ossenkopf & Henning (1994) dust opacities 
(Table 1, Column 5, commonly referred to as 0H5 dust) 
to 271.1 GHz (the effective central frequency of the Bolo- 
cam bandpass convolved with a 20 K blackbody modified 
by an opacity varying with frequency as csc i^^-^; see 
Aguirre et al. 2010) and find k,. = 0.0114 cm^ g-i. 

We make the simplifying assumption that the dust and 
gas are collisionally coupled such that Td = Tk and use 
the derived Tk in calculating Mjso. As noted in Section 
16.21 the derived Tk are likely an upper limit to the aver- 
age Tk for an entire BGPS source. Thus, the assumption 
that Td = Tk characterizes the entire BGPS source will 
overestimate Td and underestimate Miso- As discussed in 
DIO, based on NH3 maps of molecular clouds presented 
in Zinchenko et al. (1997), assuming the peak Tk for the 
entire source may underestimate Miso by a factor of up 
to 2. Additionally, upper limits due to non-detections of 
the NH3(2,2) line provide only lower limits on Miso- 

Miso could also be overestimated if some of the 1.1 
mm emission is radio free-free emission rather than ther- 
mal dust emission. Dunham et al. (2011) found that the 
mean contribution of free-free emission is only 3% based 
on a comparison of BGPS flux densities and 6 cm radio 
emission from Urquhart et al. (2009). For most BGPS 
sources, free- free emission will not significantly increase 
the flux density, but for a small number of sources ap- 
proximately half or more of the 1 . 1 mm flux density could 
be a result of free-free emission. 
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Fig. 13. — The distribution of isothermal masses Miao(120") (top 
panel) and Miao(int) (bottom panel) for the full sample (gray his- 
tograms) and the Tk subsample (striped histograms). 

Miso(120") and Miso(int) are listed in Table [71 The 
distributions of Afiso(120") and Miso(int) are shown in 
Figure [T5] for the full sample (gray histograms) and the 
Tk subsample (striped histograms). The Tk subsample 
is characterized by {log{Misoil20") /Mq)) = 2.95 ± 0.61 
and {log{Misoiint)/MQ)) = 2.94 ±0.66. See Table[5]for 
additional statistics for these mass distributions. The 
lowest mass sources (Miso< 10 Mq) at distances less than 
1 kpc, and with velocities near km have unreliable 
kinematic distances as discussed in Section [5] 

We also calculate the virial mass for each source as- 
suming a spherical, uniform density gas cloud given by 

AU,, = (5) 

where ayisr is the Gaussian line width, R is the physical 
radius, and G is the gravitational constant. This equa- 
tion assumes that the NH3 line width observed at the 
peak of the 1.1 mm emission is also valid at the edge of 
the observed 1.1 mm emission. This assumption is un- 
likely to hold for the larger BGPS sources at the farthest 
distances, as NH3 emission is not likely to trace entire 
molecular clouds but rather the denser regions within the 
clouds. As discussed in DIO, these assumptions likely re- 
sult in an overestimate of M^\^. Zinchenko et al. (1997) 
found that the line width decreases away from the peak 
of the NH3 emission in approximately half of their sam- 
ple, which would lead to an overestimate of Mvir. Ad- 
ditionally, the BGPS sources are not likely spherical or 
uniform density. The mean aspect ratio of the full sample 
is 1.5 ± 0.3. We cannot characterize the density distri- 
bution of the BGPS sources with our current data set, 
so we consider the mean power-law density profile found 
from modeling dust emission from massive star-forming 
regions, (p) = 1.8±0.4 (Mueller et al. 2002). If we correct 
for shape and non-uniform density as described in DIO, 
(Mvir/Miso) would be decreased by 36%. Since we do not 
have a measure of the density power-law, or NH3 maps 
to quantify these effects for each BGPS source, we do 
not include these corrections in subsequent analysis but 
note that the virial masses presented are upper limits. 

Mvir(int) is given in Table [T] The full sample is char- 
acterized by (log(Mvir/M0)) = 2.85 ± 0.66 while the Tk 
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Fig. 14. — Mvir/Miso(int) ratio versus Miso(int) for the full sam- 
ple. The dotted line marks the mean virial parameter of the full 
sample, (Mvir/Miso) = 1.82, and the dashed line marks the median 
virial parameter of 0.74. 

subsample is characterized by (log(Aivir/M0)) = 2.87 ± 
0.64. Figure [T4l plots the virial parameter, Mvir/Mjso, 
versus Miso(int). The full sample is characterized by 
(log(Mvir/Miso)) = -0.09 ± 0.45, and a median of 
2Q-o.i3^ which is more representative of the full sam- 
ple since the mean is skewed by the few sources with 
a virial parameter of ~100. These BGPS sources with 
Afvii/Afiso~100 are the local sources with unreliable kine- 
matic distances and upper limits for the source sizes. 
Since the virial parameter is proportional to , the ob- 
served scatter in the virial parameter could be a result of 
errors in the derived kinematic distances. Additionally, 
the scatter could be due to deviations from a uniform 
density distribution or varying source ellipticities. Alter- 
natively, the assumption that the observed line width is 
applicable at the boundary of the 1.1 mm emission could 
be invalid. 

Afiso(int) is shown versus kinematic distance in Figure 
112b . The gray circles mark the Tk subsample, and the 
gray crosses mark BGPS sources with an upper limit for 
Tk, and subsequently a lower limit in Afjso. The dashed 
black line marks the mass within a beam corresponding 
to a 5cr detection, where we assume a = 0.020 ± 0.006 
Jy beam~^ from the mean and standard deviation of 
the rms in each BGPS field, and Td= 20 K. The black 
squares and black line connecting them mark the mean 
Afiso within 2 kpc wide bins. The dotted lines mark lim- 
its between clouds, clumps, and cores of 27.5 Mq and 
750 M0 based on the following characteristic values of 
each type: clouds have M — 10'^ — 10^ Mq, clumps have 
A/ = 50 - 500 Mq, and cores have Af = 0.5 - 5 Mq 
(Bergin & Tafalla 2007). The mean Afjso increases with 
distance, and all masses are above the limiting mass per 
beam because the sources are significantly extended com- 
pared to the beam. 

6.4. Sound Speed and Non-thermal Motions 

The observed line widths provide a measure of the 
internal motions within each BGPS source. As dis- 
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T,„ [K] M,J]nt) [0] 

Fig. 15. — Left: The non-thermal Hne width versus Tk for the full sample. Black circles mark the Tk subsample and the arrows mark 
sources with only an upper limit for Tk due to a non-detection in the (2,2) line. The black line denotes the thermal sound speed given by 
a = {kTji/ ^mjf)^^^ . The non-thermal component is, on average, three times larger than the thermal sound speed. Right: Mach number 
{cviarl x/fcTk/M'") versus i\/iso(int). Here, crosses mark sources without a (2,2) detection. 



cussed in Section |4j we observe velocity dispersions of 
up to 4 km s~^, with 98% of sources having avisr < 
2.0 km s^^ , which is similar to line widths observed 
toward other samples of massive star-forming regions 
(e.g. AT/(FWHM) ~ 1 - 4 km s'^; Anglada et al. 1996) 
but an order of magnitude larger than seen in NH3 obser- 
vations of cores in nearby, low-mass star-forming regions 
(e.g. Rosolowsky et al. 2008). 

The observed velocity dispersion is a result of both 
thermal and non-thermal motions of the gas. Given the 
derived Tk we can calculate the thermal contribution 
to the line width and remove it to determine the non- 
thermal velocity dispersion: 

"^^^y^^'-^lS' 

where (fcTx/lTmi/)^/^ is the thermal broadening due to 
NH3, and niH is the mass of a single hydrogen atom. 
The full sample can be characterized by (ctnt) < 0.75 ± 
0.49 km s~^, and the Tk subsample by (ctnt) — 0.81 ± 
0.53 km s~^. The non-thermal velocity dispersions for 
the full sample are shown in Figure [15] versus Tk . The 
Tk subsample is shown as black circles, while sources 
without a (2,2) detection and an upper limit in Tk are 
plotted as arrows marking the upper- limit in Tk- The 
solid black line depicts the thermal sound speed as a 
function of Tk given by a = (fcTK//im//)^/^, where /i = 
2.37. For most sources, the non-thermal line width is 
significantly larger than the thermal component; we find 
a mean ratio of non-thermal to thermal sound speed of 
3.2 ± 1.8 for the fuh sample and 3.2 ± 1.7 for the Tk 
subsample. The sound speed and non-thermal velocity 
dispersion are listed in columns 5 and 6 of Table [51 

The right panel of Figure [121 plots the Mach number 
(given as uvisr/o) versus A'/iso(int). The non-thermal 
velocity dispersion weakly scales with mass. As discussed 
in the previous section, Mjso increases with distance, and 
the increase in ctnt seen with mass is likely the result of 



observing larger objects at greater distances. We find a 
mean Mach number of 3.23±1.87 for the full sample, and 
3.24±1.73 for the Tk subsample. 

6.5. Volume and Surface Densities 

We calculate the volume-averaged density based on 
Miso(int) and R as n = 3Miso(int)/47rR^. This quantity 
is listed in Column 6 of Table [7j and shown in Figure 
\Vlb . The dashed line marks the volume-averaged density 
corresponding to the limiting mass within a beam size 
shown in panel b and described in Section 16.31 In con- 
trast to the limiting mass, this measure of density is not 
a hard limit. It is possible to have lower and higher den- 
sity BGPS sources. For example, a BGPS source with an 
extracted size equal to the beam size will have a higher 
density if the flux is greater than the 5cr limit required 
for detection. A large, extended BGPS source with an 
extracted size much larger than the beam size can have 
a density lower than this limit. The observed flux is 
roughly linear with source size while the volume is pro- 
portional to R'^ , thereby allowing a larger source to have 
a lower density (see Section 7.2 of DIO). This trend can 
be seen in Figure [T^ where the volume-averaged density 
decreases with kinematic distance while the extracted 1 . 1 
mm source size increases. The dotted lines in this fig- 
ure denote the typical range of clouds, clumps and cores 
(clouds, n — 50 — 500 cm~'^; clumps, n = 10^ — 10^ cm~^; 
cores, n = 10'' - lO'"^ cm^^; Bergin & Tafalla 2007). 

The full sample has a minimum volume-averaged den- 
sity of 120 cm~'^, a maximum of 8.0 x 10^ cm~'^, and 
(log(n/cm-3)) = 3.12 ± 0.56. The Tk subsample has 
(log(n/cm"3)^ ^ 3 27 ± 0.58. Similarly, we calculate 
the average surface density as S = Misoiint) /irR^ . This 
quantity is listed in column 8 of Table [Tj The distribu- 
tion of surface density for the full sample ranges from 
0.006 g cm-2 to 0.24 g cm^^ with (S) = 0.037 ± 0.030 g 
cm~^. The full statistical characteristics of each sample 
are given in Table [8] 
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Fig. 16. — Observed T^y,{1,1) versus kinematic distance for 
the full sample. Circles mark the Tk subsample and crosses 
mark sources without a (2,2) detection. The decrease in observed 
TinbCljl) with distance is indicative of beam dilution. 



Fig. 17. — Excitation density versus 1.1 mm volume-averaged 
density. Black circles mark the Tk subsample while crosses mark 
sources with only an upper limit for Tk. The dashed line denotes 
a one-to-one correlation and is not a fit to the data. 



6.6. Excitation Temperature and Density 

The observed Tinb(l,l) are shown versus kinematic dis- 
tance in Figure I16[ where circles mark the Tk subsam- 
ple and crosses mark sources without a (2,2) detection. 
There is a clear decrease in Tmb(l,l) with distance sug- 
gesting that distant sources suffer from beam dilution. 
Such beam dilution causes Tex to be underestimated, 
which in turn underestimates the excitation density as 
discussed below. Such beam dilution will not affect the 
calculated Tk since it is based on the ratio of (1,1) and 
(2,2) emission, which will suffer from the same beam di- 
lution. 

The measured excitation temperatures are listed in Ta- 
ble [S] Sources not well fit by the assumed NH3 model 
have arbitararily been set to have Tex = Tk- The full 
sample is characterized by (Tex) = 4.1 ± 2.2 K, and the 
Tk subset by (Tex) = 4.5 ± 2.1 K. The mean difference 
between the excitation temperature and kinetic temper- 
ature is (Tk — Tex) = 11-5 ± 5.1 K for the full sample. 
This difference can be explained by a combination of the 
low volume-averaged densities detected and the effects 
of beam dilution. The densities are low enough that the 
NH3 transitions are not thermalized, as expected for the 
mean density we find (see DIO for more detail). 

We can calculate an excitation density, defined to be 
the density which produces the observed level popula- 
tions described by Tex for the observed Tk, from Tex, 
Tk, and r. 



where 



T„ 



fc( J(Tex) - JiT,mb)) 
/IZ.(1,1)(1- J(Tex)/J(TK)) 



J{T) 



/s(l_e-'>i'(i,i)/fcT)' 



(8) 



2.73 K, the escape probability is given by /? 



(1 — e ^)/r, where we take r = T(i^i) x 0.233, which is 
the maximum optical depth in a single hyperfine line of 



the NH3(1,1) transition. Since Tex will be lowered due to 
the beam dilution in the more distant sources while Tk is 
unaffected, the excitation density will be underestimated 
in sources suffering from beam dilution. 

The excitation densities are listed in Column 7 of Table 
[71 We cannot calculate riex for the sources where we set 
Tex = Tk. The full sample has (log(nex/cm-3)) = 3.10± 
0.35. The excitation density is plotted versus the volume- 
averaged 1.1 mm density in Figure [T71 where black circles 
mark the Tk subsample and crosses mark BGPS sources 
without a (2,2) detection. The dashed line marks a one- 
to-one correlation and is not a fit to the data. For the 
full sample, we find (log(ne2;/7i)) = —0.02 ± 0.58 and a 
median of 10*^ °^, suggesting that both density estimates 
are reasonable. 

6.7. Column Densities and NH3 Abundance 

We calculate two measures of the H2 column density, 
the column density in a beam centered on the peak of 
the 1.1 mm emission [N^^"^) and the average column 
density of the entire 1.1 mm source {N'^_^). The column 
density in a beam is calculated as 



5.(40") 



H2 



r!(40")Mff.TOHK.B,(TD)' 



(9) 



where 5*1/(40") is the flux density in an aperture of diam- 
eter 40" without the aperture correction of 1.46 required 
to correct for point source emission that falls outside of 
the aperture, f2(40") = 2.95 x 10~* sr is the sohd angle 
of the 40" aperture, ^mh-^ = 2.8, m.H is the mass of a 
single hydrogen atom, is the dust opacity described 
in section 15751 and ByiT^:,) is the Planck function evalu- 
ated at Tq. We use 5.(40") as a measure of the average 
1 . 1 mm surface brightness to measure 

j^beava rather than 

the single peak 1.1 mm pixel in order to reduce the af- 
fects of pixel-to-pixel variations caused by the noise in 
the BGPS images. The solid angle subtended by the 40" 
aperture is very close to the solid angle of the 33" effec- 
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tive beam, Q.{beam) = 2.9 x 10 * sr. We calculate N'^^^ 
as Miso(int)//i//2m//7ri?^ and note that this quantity is 
typicaUy smaher than N'^^™- since it is averaged over the 
entire area ol the BGPS source while N'^^"' is calculated 
at the peak of the 1.1 mm emission. M^e^™ and N"^^^ 
are listed in columns 9 and 10 of Table [T] 

The average H2 column densities, N'j^^ , are shown in 
Figure versus kinematic distance. The gray circles 
denote the Tk subset, gray crosses denote sources with- 
out a (2,2) detection, and the black squares and black 
solid line denote the mean H2 column densities in 2 kpc 
bins. The dashed line denotes the H2 column density 
corresponding to the mass per beam limit described in 
section [^751 The dash-dot line marks a threshold gas sur- 
face density of 122.5 M0 pc~^ above which efficient star 
formation occurs (116 M© pc-^ Lada et al. 2010; 129 
M0 pc~^, Heiderman et al. 2010), and the dotted line de- 
notes a surface density of 1 g cm~^ required to prevent 
excessive fragmentation and allow massive stars to form 
(Krumholz & McKee 2008). See section for further 
discussion of these limiting column densities. The mean 
H2 column density decreases slightly with kinematic dis- 
tance, but is essentially flat. The full sample is described 
by (log(iV^7/cm-2)) = 21.81 ±0.27, and 180 M© pc-^. 

Using Rv = 3.1 and A^^^/Ay = 9.4 x lO^" cm'^ mag-i 



(Bohlin et al. 1978), the mean N'^^ for the full sample 
corresponds to 8.5±6.8 mag of visual extinction. Assum- 
ing a Rv = 5.5 extinction model (Weingartner & Draine 
2001), Nl^iyAv = 6.86 x lO^o cm'^ mag-^, the mean 
of the full sample corresponds to 11.7 ± 9.3 mag. 

We calculate the NH3 column densities as described in 
DIG and list them in column 11 of Table [T] Statistical 
properties are listed in Tabled 

The NH3 column densities are described 
by (log(iV7Vff3/cm-2)) = 14.75 ± 0.37 and 
(log(iVjvff3/cm-2)) = 14.95±0.29 for the fuh sample and 
Tk subsample, respectively. We additionally find a mean 
NH3 abundance of (log(iVjvff3/iV|f7")) ='-7.34 ±0.03 
for the full sample. Our measured NH3 abundance 
is consistent with other studies: 3 x 10~* (Harju et 
al. 1993), 2.8 x 10"^ (Tafalla et al. 2006), and 2 x 10^^ 
(Foster et al. 2009). 

7. GALACTIC TRENDS IN PHYSICAL PROPERTIES 

As discussed in section [5l each BGPS source has a 
unique Rcai given by the observed radial velocity and 
assumed Galactic rotation curve, allowing us to study 
the trends in physical properties as a function of distance 
from the Galactic center. In this section we include the 
Gemini OBI molecular cloud (Gem OBI) sources from 
DIO as an outer Galaxy data point. Since the BGPS 
did not perform a blind survey of the outer Galaxy, but 
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Fig. 19. — Mean properties of near BGPS sources versus Galactocentric radius for (a) kinematic distance, (b) Tk, (c) Afiso, (d) i?, (e) n, 
(f) N"^^ , {g)N{^[f^, and (h) Nj^u^/N^^"^ . Each point represents the mean over 39 BGPS sources, except for the last two points, which 
represent 38 and 33 sources, respectively. The points are centered at the mean Rq^i of 3.3, 4.2, 4.6, 4.9, 5.5, 7.0, and 10.5 kpc. The blue 
diamonds and vertical lines denote the mean and standard error in the mean (given by a/y/N where a is the standard deviation and TV 
is the number of sources per bin) for each property. The blue horizontal lines in panel (a) represent the range of Rq^i included in each 
point. The dotted line in panel (a) represents the distance to the nearest source at a given -Rcali which correspnds to the 1 = 0° line of 
sight. The dashed lines represent: (c) the 50(t mass per beam, (d) beam size, which represents a lower limit on extracted source size, and 
the scale above which we filter out uniform emission representing an upper limit to extracted source size (5.9', dash-dot line), (e) the 5a- 
volume-averaged density per beam, and (g) the 50o" average column density per beam. 



rather targeted known star-forming regions (Aguirre et 
al. 2011), the Gem OBI sources may be biased in terms 
of evolutionary stage. We discuss possible effects of this 
bias in Section 17.21 We also note that the Gem OBI 
sources may not be representative of all outer Galaxy 
star formation. 

Figure US] plots the physical properties as a function 
of Galactocentric radius for the Tk subsample (gray cir- 
cles) and the BGPS sources without a (2,2) detection 
(gray crosses) . The black squares denote the mean prop- 
erties within 2 kpc bins centered at 1, 3, 5, 7, 9, and 
11 kpc. We observe BGPS sources at both the near and 
far kinematic distances for each Galactocentric radius. 
The concentration of BGPS sources in lines in Figure 
[TSk is due to the range of Galactic longitudes included 
in the full sample. Each sideways V formation corre- 
sponds to one of the four ranges in Galactic longitude, 
where the vertex of the V created by the ^ ~ 9° sources 
is at the smallest i?Gai and the I ^ 54° sources are found 
at i?Gai~ 7 kpc. 

The observed dichotomy of kinematic distance as a 



function of Galactocentric radius causes a spread in phys- 
ical properties such as Mjso, n, and as a function of 
RgbX because of their dependence on the kinematic dis- 
tance (see Figurc[T2|l. This spread in properties will serve 
to mask any trends in properties as a funciton of i?Gai- In 
this section, we consider the near subsample (233 BGPS 
sources placed at the near kinematic distance) in order 
to remove the ambiguity introduced by the kinematic 
distance dichotomy. 

The mean and standard error in the mean are pre- 
sented for each property for the near subsample in Figure 
[TOl The 233 BGPS sources were sorted by i?Gai and spHt 
evenly into six bins with 39 sources in the first five bins 
and 38 sources in the sixth bin. The 33 Gem OBI sources 
are presented separately in the seventh bin. This binning 
method provides more measurements of the mean prop- 
erties where there are more BGPS sources. The 7 bins 
are centered at i?Gai of 3.3, 4.2, 4.6, 4.9, 5.5, 7.0, and 
10.5 kpc. Bins 2-4 centered at 4.2, 4.6, and 4.9 kpc 
represent sources in the 5 kpc molecular ring. The blue 
diamonds represent the mean properties within each bin 
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and the vertical blue lines represent the standard error in 
the mean, which is given by cr/V7V where a is the stan- 
dard deviation and N is the number of sources in each 
bin. 

Figure [T9k plots the mean kinematic distance versus 
i?Gai ■ The horizontal blue lines depict the range of -Rcai 
values included in each bin. The dotted lines correspond 
to the distance to the nearest source at a given -Rcai, 
which lies along ^ = 0° between the sun and the Galactic 
center. The near subsample decreases in kinematic dis- 
tance as i?Gai increases such that sources at higher i?Gai 
are at smaller kinematic distances. The kinematic dis- 
tances are similar to the nearest source distance shown 
by the dotted lines allowing us to apply our understand- 
ing of the biases as described in Sections [6l and [8T2] to the 
near subsample as a function of i?Gai- 

7.1. Observed Trends 

Here we highlight the observed trends in physical prop- 
erties with i?Gai, and discuss the possible biases in the 
following section. 

The sources within the molecular ring are located 
at larger kinematic distances than expected given the 
smooth trend of decreasing kinematic distance with in- 
creasing i?Gai seen in the inner Galaxy sources located 
outside of the molecular ring (see Figure [T9k ) . We will 
explore the possible effects of the larger kinematic dis- 
tance toward the molecular ring sources in the following 
section. 

Figure [T9b shows the kinetic gas temperature derived 
from the NH3 observations. There is no real trend in 
Tk of sources in the inner Galaxy, but there is a 3 K 
difference between the mean Tk of the inner Galaxy and 
Gem OBI sources. 

Afiso is shown in Figure [T9b . The dashed lines depict 
a 50(7 mass per beam at distances corresponding to the 
dotted lines in panel (a). Mjso steadily decreases with 
Roai with the trend roughly following the shape of the 
curve representing a 50cr mass per beam. The three bins 
closest to the molecular ring at 4.2, 4.6, and 4.9 kpc 
deviate from the general trend in decreasing Mjso with 
i?Gai and are roughly a factor of a few more massive than 
described by the general decreasing trend. 

Figure [T9H displays the mean observed physical radius 
in each bin. The dashed lines mark the beam size at 
distances corresponding to the dotted lines in panel (a), 
which represent a minimum source size, and the dash-dot 
lines mark the upper limit in source size above which the 
data reduction filters out emission. The mean radius de- 
creases with i?Gai , but two of the three bins representing 
the molecular ring have radii approximately 1.5 times 
that expected from the general decreasing trend. 

The volume-averaged density is shown in Figure [T9b as 
a function of -Rcai- The dashed lines mark the 5cr density 
per beam required for a detection in the BGPS. Two of 
the bins representing the molecular ring fall below this 
line while the other points lie above. 

Figure [TW plots N^^^ as a function of i?Gai • The av- 
erage H2 column density is roughly constant with i?Gai, 
with the lowest column densities seen in the molecular 
ring and the highest in Gem OBI. The dashed line de- 
notes a 50cr column density per beam. The observed NH3 
column density is shown in Figure 119b as a function of 



^Gai- The column density decreases with i?Gai and only 
one of the molecular ring bins falls below the general 
trend. Figure [T9h plots the observed NH3 abundance 
versus i?Gai- The abundance, given by NMH^/N'jf^™, 
decreases with i?Gai and is approximately a factor of 7 
higher in the inner Galaxy than in the Gem OBI sources. 

7.2. Biases and Implications 

As mentioned above, the Gem OBI sources belong to 
a molecular cloud known to be forming stars. This could 
bias the physical properties thereby providing an inap- 
propriate comparison point for the inner Galaxy sources. 
The Gem OBI sources differ dramatically from the in- 
ner Galaxy sources in only two properties: Tk and NH3 
abundance. The warmer mean Tk observed in the Gem 
OBI sources could be due to a bias toward later evo- 
lutionary stages, which typically have warmer dust tem- 
peratures (Battersbyet al. 2010). Similarly, Tk of BGPS 
sources increases with the number of associated mid-IR 
sources (see section [873]) . By targeting a region of known 
star formation, we have potentially excluded the earli- 
est evolutionary stages that have colder temperatures, 
therefore resulting in a warmer mean Tk. 

Foster et al. (2009) found no significant difference be- 
tween the NH3 abundance of starless and protostellar 
cores in Perseus, suggesting that the potential bias of 
the Gem OBI sources to later evolutionary states is not 
likely to be the cause of the lower NH3 abundance in 
Gem OBI sources. The general decrease in NH3 abun- 
dance (Figure [TSh) could be the result of changes in dust 
properties and/or mctallicity with i?Gai- Changes in the 
gas-to-dust ratio, grain size distribution, and tempera- 
ture of dust grains across the Galaxy could all result 
in the observed decrease in NH3 abundance with i?Gai- 
Similarly, a decrease with i?Gai in the amount of nitro- 
gen available to form NH3 would result in a lower abun- 
dance at larger i?Gai- The decrease we observe mimics 
the decreasing nitrogen abundance seen as a function of 
i?Gai (e.g. Gummersbach et al. 1998; Shaver et al. 1983; 
Smartt et al. 2001; RoUeston et al. 2000). We fit a line 
to the NH3 abundance as a function of i?Gai and find 

JNI13 ■ ^ 

We observe a gradient of -0.096 dex/kpc in the NH3 
abundance, which agrees within the error bars with 
the observed nitrogen gradient (—0.09 ± 0.015 dex/kpc; 
Shaver et al. 1983). Thus, the NH3 abundance decreases 
at least as quickly as the nitrogen abundance but possibly 
faster due to possible changes in dust properties to which 
our observations are not sensitive. This trend needs to be 
taken into account when calculating H2 column densities 
and masses based on observed NH3 column densities. 

The effects of distance must also be taken into account 
when intrepreting trends with i?Gai • The Malmquist bias 
together with the spatial filtering in the BGPS data has 
resulted in the detection of different types of structures 
at different distances (i.e. cores, clumps, clouds; see sec- 
tion [821), which complicates the study of BGPS sources 
as a function of position in the Galaxy. This complica- 
tion, which will affect all continuum studies, needs to be 
addressed carefully. 

Both the mass and radius of BGPS sources exhibit the 
same general trend of decreasing values with increasing 
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Fig. 20. — Distributions of kinematic distance (upper panel), physical radius (lower left panel), and A/iso(int) (lower right panel) for 
the near kinematic distance sources located within the molecular ring (solid line), outside the molecular ring (light gray), and outside the 
molecular ring with 3.5 kpc< D < 6.5 kpc (dark gray). The legend gives the KS probabilities comparing the full sample of sources inside 
and outside of the molecular ring (light gray) and for the kinematic distance selected samples of sources inside and outside of the molecular 
ring (dark gray). 



^Gai with peaks within the molecular ring. For the near 
subsample, i?Gai and D are inversely proportional so that 
large values of i?Gai correspond to small kinematic dis- 
tances and small i?Gai corresponds to large kinematic 
distances. Therefore, the larger masses and radii seen 
at low i?Gai are likely the result of the Malmquist bias. 
Indeed, the slope of the mean mass versus i?Gai roughly 
follows the slope of the 50cr mass per beam line suggest- 
ing that trend is caused by distance effects. 

The general trends of decreasing masses and radii 
with increasing i?Gai are due to distance effects, but are 
the larger masses and radii seen toward molecular ring 
sources real? The difference in mean masses and radii 
plotted in Figure [12] suggest the molecular ring sources 
are truly larger and more massive, but the sources with 
-RGai= 4 — 5 kpc are also located at larger distances 
than described by the general trend of decreasing kine- 
matic distance with i?Gai (Figure [19^). We explore the 
possible effects of distance by comparing properties of 
sources inside (4 < i?Gai < 5 kpc; 115 BGPS sources) 
and outside the molecular ring (151 BGPS sources) using 
Kolomogorov-Smirnov (KS) tests. The two-sided KS test 
compares the cumulative distribution functions of two 
samples and measures the probability that the sources 



are drawn from the same parent distribution, but high 
KS probabilities do not necessarily indicate the same par- 
ent population. Figure!^ plots the distributions of kine- 
matic distance (upper panel), physical radius (lower left 
panel), and Miso (lower right panel) for sources within 
the molecular ring (solid line histogram) and sources out- 
side the molecular ring (light gray histogram). The KS 
probability for the kinematic distance (10~^^) shows that 
we are not comparing similar distance distributions by 
splitting the sample into sources found inside the molec- 
ular ring and those found outside the molecular ring. The 
sources located outside the molecular ring span a large 
range of kinematic distances, as seen in the upper panel 
of Figure [20l We have attempted to further untangle 
the effects of distance by considering only the sources 
inside and outside the molecular ring within the same 
range of kinematic distances, specifically 3.5 < D < 6.5 
kpc. This distance range includes 114 of the 115 sources 
inside the molecular ring, and 42 of the 151 sources out- 
side the molecular ring (dark gray histogram in Figure 
[20]) . We performed similar KS tests on these kinematic 
distance selected samples and found KS probabilities of 
0.41, 0.081 and 0.036 for D, R, and Mi^o, respectively. 
The KS probabilities have increased compared to the 
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non-kinematic distance selected sample and no longer 
suggest measurable differences for any of these proper- 
ties. 

Therefore, the larger mean masses and radii of sources 
within the molecular ring are likely caused by distance 
effects. When we compare sources at similar kinematic 
distances, we see that there is little variation in physical 
properties regardless of environment (inside versus out- 
side the molecular ring, for example). This is in contrast 
to the molecular clouds seen in the ^^CO GRS survey. 
The molecular clouds within the molecular ring typically 
have warmer temperatures, larger areas, larger masses, 
higher peak H2 column densities, and typically contain 
more clumps than molecular clouds located outside the 
molecular ring (Rathborne et al. 2009). This suggests 
that environment may affect the large scale molecular 
clouds while the regions of dense gas (n ^ lO'^ cm~'^) 
within the molecular clouds exhibit similar properties re- 
gardless of environment. 

8. DISCUSSION 

8.1. Comparison to Other Surveys 

Schlingman et al. (2011) conducted a survey of BGPS 
sources in HCO+(3-2) and N2H+(3-2) and were able to 
break the kinematic distance ambiguity for 648 sources 
that they placed at the near distance. They find a me- 
dian radius of 0.75 pc, median isothermal mass assuming 
Td= 20 K of 330 M0, and a median volume- averaged 
density of 2.5 x 10^ cm'^. For the 233 BGPS sources 
placed at the near kinematic distance in our sample, we 
find a median radius of 0.95 pc, median mass of 350 
Mq, and a median volume-averaged density of 2.2 x 10'^ 
cm~^, very similar to the properties of the Schlingman 
et al. sample. In contrast to our sample, where we find 
a median virial parameter of 0.74, Schlingman et al. find 
typical virial parameters of a few and only a small num- 
ber of sources have virial parameters less than one. This 
difference is likely due to the larger line widths detected 
in HCO+ (median FWHM= 2.98 km s"^) compared to 
the line widths detected in NH3 (median FWHM = 0.64 

km s"^ %/8h^ = 1.5 km s'^). 

The APEX Telescope Large Area Survey of the Galaxy 
(ATLASGAL; Schuller et al. 2009) has surveyed the 
Galactic plane at 870 ^m. Early results from this survey 
and NH3 follow-up observations reveal similar physical 
properties as those presented in this work. For 24 AT- 
LASGAL sources near £ ~ 19° that were detected in 
NH3, Schuller et al. find Tk ranging from 14 to 33 K, 
masses ranging from roughly 10^ to 10'* Mq, and H2 col- 
umn densities from 7x 10^^ cm^^ to 1.3x 10^^ cm~^. The 
range of properties of ATLASGAL sources are consistent 
with the range of properties of BGPS sources presented 
here. 

Compared to the GRS clouds studied by Roman-Duval 
et al. (2009), the BGPS sources are smaller (BGPS R = 
0.01- 8.5 pc; GRS i? = 0.1-40 pc), less massive (BGPS 
Miso= 0.1-2.4x10^ Mq; GRS M = 10-10^ Mq), higher 
density (BGPS {n)= 5x10^ cm^^. Qj^g (^^^^^ 230 cm'^), 
and have higher surface densities (BGPS (S)= 180 Mq 
pc~2; GRS (1;)= 144 Mq pc"^). This is consistent with 
the paradigm that the BGPS sources are typically the 
smaller, denser regions located within molecular clouds. 




Distance 



Fig. 21. — Fraction of sources in the full sample of 456 BGPS 
sources detected in NH3 nearer that are cores (thin black lines) and 
farther that are clouds (thick blue lines) as a function of kinematic 
distance determined by radius (solid lines), mass (dotted lines), 
and volume-averaged density (dashed lines). Both ends of the d is- 
tance range suffer from small number statistics. See section [52] for 
further details. 



8.2. What IS a BGPS Source? 

DIO addressed the question "what is a BGPS source" 
in terms of expected masses and densities as a function 
of distance based on sources at a single, well-known dis- 
tance. With a large sample spanning a large range of 
kinematic distances, we can address what type of source 
the BGPS is detecting at each distance empirically. As 
discussed throughout section \6\ BGPS sources can be 
classified as clouds, clumps, or cores based on their size, 
mass, and volume-averaged density. While there can be 
significant overlap in the physical properties of clouds, 
clumps, and cores, they can be approximately separated 
by the dotted lines shown in Figure ll^ a-c'). 

Figure [5T] shows the fraction of sources closer than a 
given distance that are described as cores (thin black 
lines; R < 0.125 pc, Mi^o < 27.5 Mq, and n > 10'^ cm'^) 
based on radius (solid), mass (dotted), and volume- 
averaged density (dashed). We find that greater than 
90% of the BGPS sources out to a distance of 0.68 kpc are 
well-described as cores. However, we suffer from small 
number statistics at distances less than 2 kpc (see Figure 
[9|), and the distance out to which we detect cores could 
be greater. Since no sources are detected between 1 and 
approximately 1.75 kpc, the fraction of BGPS sources 
that are cores can not be characterized in this distance 
range. Thus, this result is consistent with DIO, which 
determined that cores are detected out to approximately 
2 kpc. 

Similarly, we determine the fraction of sources lo- 
cated farther than a given distance that are character- 
ized as clouds {R > 1.25 pc, Miso > 750 Mq, and 
n < 750 cm~^). These fractions are shown in Fig- 
ure [21] as the thick blue lines. The fraction of sources 
with volume-averaged densities below the limit for clouds 
never reaches 0.9, but does follow the general trends seen 
in the fractions based on radius and mass. The density 
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Fig. 22. — Fraction of BGPS sources associated with at least one mid-IR source (Dunham et al. 2011) as a function of kinematic distance 
(left) and iJcal (right) for the near subsample of 233 BGPS sources placed at the near kinematic distance. Each diamond represents 
includes 39 BGPS sources, except for the bin at Rq^i^I kpc, which represents 38 BGPS sources. 



curve never reaches 90 percent because it is suffering from 
effects of volume-averaging. Thus, we do not consider the 
volume-averaged density when determining the distance 
at which we detect >90% clouds. We find 7.3 kpc based 
on the fraction with radii > 1.25 kpc, and 6.6 kpc from 
the fraction with Miso(int) > 750 Mq. We take the aver- 
age of these two distances, 7.0 kpc, as the distance above 
which the BGPS detects >90% clouds. 

At distances up to approximately 1 kpc, the BGPS 
detects cores. At distances greater than approximately 
7 kpc, the BGPS detects clouds. At intermediate dis- 
tances, the BGPS will detect predominately clumps, but 
also cores and clouds. While the high column density fea- 
tures seen in the BGPS do identify regions of dense gas, 
they represent very different scales of structures when 
looking through the Galaxy. Similar continuum based 
studies (e.g. ATLASGAL, Schuller et al. 2009; Hi-Gal, 
Molinari et al. 2010) will also be subject to these effects. 
Kinematic data is crucial in interpreting the results of 
such studies and determining what structures are actu- 
ally observed. 

8.3. Star Formation Activity 

Dunham et al. (2011) characterized the star forma- 
tion activity of the BGPS sources based on the presence 
of mid-infrared (mid-IR) sources along a line of sight 
coincident with each BGPS source. They placed each 
BGPS source in one of four groups representing increas- 
ing confidence in the type of associated mid-IR source 
reliably indicating star formation activity (as opposed 
to a chance alignment). Group 3 represents the highest 
confidence in star formation activity and includes BGPS 
sources matched with an Extended Green Object (EGO; 
Cyganowski et al. 2008) or Red MSX Survey (RMS; 
Hoare et al. 2004; Urquhart et al. 2008) source. Groups 
1 and 2 include BGPS sources matched with GLIMPSE 
red sources cataloged by Robitaille et al. (2008) or a 
deeper GLIMPSE red source catalog created by Dun- 
ham et al. (2011). Group includes BGPS sources that 
were not matched with any mid-IR sources, and we refer 
to them as "starless" although they are not necessarily 



truly starless (see Dunham et al. 2011 for more details). 

Dunham et al. (2011) find that 49% of BGPS sources 
in the inner Galaxy where all surveys overlap contain at 
least one mid-IR source. We find that 67% of the full 
sample contains at least one mid-IR source, while 33% 
are "starless." 68% of the BGPS sources without a de- 
tection in the NIl3(l,l) line were "starless," and groups 
1 through 3 comprise roughly 21%, 9%, and 3%, respec- 
tively. The BGPS sources without NH3 or any corre- 
sponding mid-IR sources may be false detections in the 
BGPS catalog. 

Figure [52] plots the fraction of BGPS sources in the 
near subsample that contain at least one mid-IR source 
as a function of kinematic distance (left panel). The 
first five data points represent 39 BGPS sources, and 
the sixth represents 38. By placing the same number of 
BGPS sources in each bin, we obtain more data points 
at distances where more BGPS sources are located. The 
fraction peaks near 85% at a distance of approximately 
3.0 kpc, and reaches a minimum of 55% at 6 kpc. The 
decrease in fraction of BGPS sources with mid-IR sources 
is likely a result of the increasing distances. At the far- 
thest distances shown here, the BGPS is dectecting the 
large scale objects, clumps or clouds, and the GLIMPSE 
and RMS surveys would only detect the most luminuous 
YSOs, which become increasingly rare with increasing lu- 
minosity and mass. Thus, the decreasing fraction we see 
with distance may simply be a reflection of a decreas- 
ing number of detected YSOs in the mid-IR catalogs. 
The fraction of BGPS sources with mid-IR sources also 
reaches a minimum at kinematic distances of roughly 4.5 
kpc. Sources in this kinematic distance bin correspond to 
sources within the 5 kpc molecular ring, and we further 
discuss the drop in fraction with star formation activity 
in the following section. 

The additional gas properties provided by the NH3 sur- 
vey presented in this work allow for further exploration 
of physical properties as a function of the star formation 
activity groups defined by Dunham et al. (2011). We 
find that the mean Tk, (Jvisr, and A^jv^jj^ increase with 
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Fig. 23. — Tk versus number of mid-IR sources associated with 
each BGPS source from Dunham et al. (2011). Plus signs are group 
sources, asterisks are group 1 sources, open squares are group 2, 
diamonds are group 3, and the large, filled squares mark the mean 
Tk for each number of mid-IR sources. The (Tk) increases with 
number of mid-IR sources up to 6 mid-IR sources, above which we 
suffer from small number statistics. 



star formation activity group such that BGPS sources 
associated with an EGO or RMS source have higher gas 
kinetic temperatures, larger observed hne widths, and 
higher NH3 column densities. The mean Tk increases 
from 13.9 K in group to 22.7 K for group 3, while 
(o'visr) increases from 0.70 km s""'^ to 1.10 km s^"'^ from 
group to 3. Similarly, (N-^-^^) increases from 5.6 x 10^** 

cm~^ for group to 1.6 x 10^^ cm~^. See Table |9] for 
further statistical characterization of each property dis- 
tribution. 

We see a trend of increasing (Tk) with increasing num- 
ber of associated mid-IR sources such that BGPS sources 
with a larger number of mid-IR sources are typically 
warmer than BGPS sources with fewer mid-IR sources 
(see Figure [531) ■ There is considerable spread in Tk with 
number of mid-IR sources, and the observed trend is a 
result of increase in the minimum Tk detected with num- 
ber of mid-IR sources. The warmest temperatures in our 
sample are seen in BGPS sources across the entire range 
of number of mid-IR sources. 

We also explore the fraction of BGPS sources with star 
formation activity as function of i?Gai- The right panel 
of Figure [22] plots the fraction of BGPS sources with 
mid-IR signs of star formation versus Galactocentric ra- 
dius for the 233 sources at the near kinematic distances. 
The binning is the same as in Figure [TOj excluding the 
sources in Gemini OBI since the GLIMPSE catalogs do 
not cover the outer Galaxy. The fraction with mid-IR 
sources is lowest at the smallest i?Gai and highest at the 
largest Rcai and exhibits a 10% decrease towards the 5 
kpc molecular ring. 

The decrease in the fraction of BGPS sources associ- 
ated with a mid-IR source seen in the molecular ring 
(Figure [221) is not indicative of less star formation ac- 
tivity within the ring. We see a clear peak in sources 
toward Rcai = 4 — 5 kpc (Figure [lU in sources that are 



Fig. 24. — Distributions of BGPS sources with (solid line) and 
without (shaded gray) star formation activity as described in Dun- 
ham et al. (2011). 



associated as well as those that are not associated with 
a mid-IR source. Although the gas and dust needed to 
form stars is more abundant in the molecular ring, the 
fraction of BGPS sources with mid-IR sources is lower. 
The decrease in the fraction of BGPS sources with mid- 
IR sources could be a result of a few different effects. The 
first possibility is the increased mean kinematic distances 
to the sources within the molecular ring. The limiting lu- 
minosity (and mass) of the mid-IR catalogs considered by 
Dunham et al. (2011) increases with distance such that 
only the most luminuous (highest mass) YSOs will be 
detected towards the most distant BGPS sources. The 
overall increase in the fraction of BGPS sources with a 
mid-IR source as a function of -Rcai is therefore a result 
of the decreasing kinematic distance with increasing Rq^i 
and the rarity of high-mass stars and YSOs. Similarly, 
since the BGPS sources within the molecular ring are 
located at larger kinematic distances, the lower fraction 
of mid-IR sources could simply be a result of the higher 
mass detection limit of the mid-IR catalogs at the larger 
kinematic distances. 

A second possible cause of the lower fraction of sources 
with mid-IR counterparts is the abundant dense gas 
within the molecular ring (e.g. Scoville & Solomon 
1975; Burton et al. 1975; Clemens et al. 1988; Kol- 
pak et al. 2002; Rathborne et al. 2009; Roman-Duval 
et al. 2010), which would provide higher extinction to- 
ward the mid-IR sources in this region. Greater ex- 
tinction would effectively raise the mass detection limit 
of the YSOs at similar distances. This could reduce 
the fraction of BGPS sources with a mid-IR counter- 
part simply because the most massive YSOs are more 
rare. The third possibility is the increased background 
of PAH emission toward the longitudes corresponding to 
the molecular ring. Point sources are easily lost in the 
increased background 8^m emission, and are therefore 
unlikely to be included in the point source catalogs. Ro- 
bitaille et al. (2008) employed a flux cut to avoid this 
issue and ensure the reliability and completeness of the 
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Fig. 25. — Miso(int) versus radius for the full sample of 456 BGPS 
sources. Symbols and colors denote the star formation activity 
within each BGPS source as defined by Dll. Black squares rep- 
resent group sources (no corresponding mid-IR sources), blue 
crosses mark group 1 sources, green diamonds denote group 2 
sources, and red circles depict group 3 sources (young RMS sources 
and EGOs). Error bars are shown for every tenth BGPS source. 
The dash-dotted line marks the Krumholz & McKee (2008) sur- 
face density threshold, the dotted line marks the Kauffmann & 
Pillai (2010) relation, and the dashed line denotes the Heiderman 
et al. (2010) and Lada et al. (2010) threshold for efficient star for- 
mation. 

GLIMPSE Red Source Catalog. However, the mid-IR 
sources considered by Dll are dominated by an addi- 
tional GLIMPSE catalog created expressly for that work, 
which did not include a flux cut and will certainly be sen- 
sitive to the background emission. It is likely that the 
reduced fraction of BGPS sources that coincide with a 
mid-IR source is the combined result of all three of the 
effects discussed here. 

8.4. Comparisons to Criteria for Efficient or Massive 
Star Formation 

In this section, we compare various criteria for effi- 
cient or massive star formation with our BGPS sample. 
We plot mass (Miso(int), Equation S]) versus radius {R, 
Equation[T]) for the full sample of BGPS sources with dis- 
tance estimates in Figure [55] Error bars are plotted for 
every tenth data point to avoid clutter. Points are coded 
by shape a nd c olor according to the four groups discussed 
in Section 18.31 Note that uncertainties in distance will 
tend to create a correlation of the form M oc . 

Heiderman et al. (2010) compared star formation rate 
surface surface densities and gas surface densi- 

ties (I](gas)) for 20 nearby clouds that are forming pri- 
marily low-mass stars and a group of more distant dense 
clumps that are forming massive stars. They found that 
T,(SF) was a strong non-linear function of T,(gas) up 
to about 129 ± 14 Mq pc~^ Above that "threshold", 
Yi{SF) was roughly linear with T,(gas). Independently, 



Lada et al. (2010) found that the star formation rate in 
a sample of clouds was linearly proportional to the cloud 
mass above Tj(gas) ~ 116 Mq pc~^, while comparisons 
to total cloud mass showed no regularity. 

In Figure [251 a dotted line of T,{gas) = 122.5M0pc-^ 
or N^^^ > 6.49 x 10^^ cm~^, represents an average of 
the criteria found by Heiderman et al. (2010) and Lada 
et al. (2010). We refer to this as the HL criterion for "effi- 
cient" star formation. Of the 456 sources in the full sam- 
ple, 211 (46.3%) satisfy the HL criterion when we use the 
surface density of the entire source. However, since the 
BGPS sources have higher surface density structures on 
smaller scales, BGPS sources that do not satisfy the HL 
criterion overall may still have subregions that do satisfy 
it. Thus, 46.3% should be viewed as a lower limit. If, for 
example, we consider the peak column density, N^°-"\ 
86.8% of the full sample of BGPS sources would have 
some part that satisfies the HL criterion. Of the BGPS 
sources with the most certain star formation (Group 3), 
70.5% (36 of 51) have average surface densities above 
the HL criterion, and 98% (50 of 51) have peak surface 
densities above the HL criterion. 

A dash-dotted line in Figure [25] shows the criterion to 
avoid excessive fragmentation of Y,(gas) = 1 g cm^^ = 
4787 Mq pc"^ advocated by Krumholz & McKee (2008), 
and presented earlier by McKee & Tan (2003). None 
of the BGPS sources in our full sample come close to 
satisfying this KM criterion. As with the HL criterion, 
parts of these sources may satisfy this criterion. The 
surveys by Plume et al. (1992) toward regions known 
to be forming massive stars found quite a few sources 
satisfying this criterion, though Wu et al. (2010) found 
that most exceeded the KM criterion only as measured in 
the highest density tracers. Since the Plume et al. sources 
are certainly included in our BGPS survey, they must 
form the densest parts of the BGPS sources. 

Kauffmann et al. (2010) and Kauffmann & Pillai 
(2010) have recently suggested a much less stringent cri- 
terion for massive star formation than the KM crite- 
rion discussed above. They argued that clouds known 
to be forming massive (M* ~ 10 Mq) stars have struc- 
tural properties described by m(r) > 870 Mq {r/pc}^-^^, 
where m(r) is the mass within radius r about a peak 
surface density. Clouds below this criterion are generally 
not forming massive stars. This "KP" criterion corre- 
sponds to a line of different slope in Figure [25] and a 
criterion for surface density that decreases with radius: 
T^{gas) (X r""-^''. While we have not characterized the 
internal structure of our BGPS sources, we can compare 
their total mass and radius to the KP relation. When 
determining their relationship, KP reduced Ossenkopf & 
Henning (1994) dust opacities by a factor of 1.5. Since 
we have assumed Ossenkopf & Henning dust opacities 
without the factor of 1.5 reduction, we scale the mass- 
radius relation to our assumed dust opacities and use 
m(r) > 580 Mq (r/pc)^-^^ for the KP relation appro- 
priate for our mass estimates. In the full sample of 
BGPS sources, 47.6% have masses above the KP rela- 
tion. Among the Group 3 sources, those most certain to 
be forming stars, and likely to be forming massive stars, 
80% (41 of 51) he above the KP relation. 

Since the BGPS includes sources ranging from cores 
to clouds, we consider all three star formation criteria 
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discussed above in order to identify all BGPS sources 
that have high enough densities to be forming stars on 
vastly different size scales. We estimate that about half 
the full sample lies above both the HL and KP crite- 
ria and should be forming stars, including massive stars, 
with high efficiency. Smaller parts of a larger fraction 
(86.8%) of BGPS sources may be forming stars efficiently, 
and a still smaller part may satisfy the KM criterion 
for massive star formation. As discussed in section 15731 
we observe mid-IR sources toward 67% of the full sam- 
ple. Since the mid-IR catalogs considered by Dunham 
et al. (2011) include high-mass YSOs across the Galaxy, 
we are capable of seeing all massive YSO counterparts 
of BGPS sources. Conversely, the mid-IR catalogs used 
in Dunham et al. (2011) include only the nearby low- 
mass YSOs. Therefore, we should not expect to detect 
mid-IR counterparts toward the full 86.8% BGPS sources 
with some surface densities high enough for efficient star 
formation. Roughly speaking, the statistics of mid-IR 
source detection are consistent with both the HL and 
KP relations. 

The 47.6% of BGPS sources that satisfy the KP rela- 
tion contain 79.3% of the total mass in the full sample. 
Thus the majority of the mass in the BGPS full sam- 
ple would be predicted by the KP relation to be form- 
ing massive stars. These sources are characterised by 
the following mean properties: (log(Miso(12O")/M0)) = 

3.29 ±0.42, (log(Miso(mt)/MQ))= 3.30 ± 0.49, {R) = 

2.30 ± 1.51 pc, (log(n(120")/cm-3))= 3.29 ± 0.49, 
(log(n(mt)/cm-3))= 3.16±0.56, (E(mt))= 0.051±0.036 
cm-2, and {\og{N^if^"^ / cm-^)) = 22.23±0.28, where the 
standard deviation of the distribution is given as an er- 
ror representing the spread in each property. The BGPS 
sources satisfying the KP relation are typically larger, 
more massive, and less dense than the previous samples 
of massive star-forming regions chosen based on tracers 
of star formation that were discussed in section [1] 

9. SUMMARY 

We present the results of a survey of NIl3(l,l), (2,2), 
and (3,3) toward a sample of 631 1.1 mm contin- 
uum sources from the Bolocam Galactic Plane Survey 
(BGPS). We have detected the NH3(1,1) line towards 
456 BGPS sources (72%), demonstrating that the high 
column density features identified in the BGPS and other 
continuum studies are excellent predictors of the presence 
of dense gas. We have determined kinematic distances 
and resolved the kinematic distance ambiguity (KDA) to 
all 456 BGPS sources via association with an 8 /im IRDC 
seen in the GLIMPSE mosaics or the presence of HIS A in 
the VGPS or SGPS. We have placed 233 BGPS sources 
at the near distance, 221 at the far distance, and two 
at the tangent distance. The kinematic distances range 
from 0.1 to 15.0 kpc. 

With the KDA resolved, we calculate physical prop- 
erties based on both the 1.1 mm continuum emission 
(e.g. R, Miso, n, E, and NH2) and the NH3 line emis- 
sion (e.g. Tk, Tex, Mvir, avisr, (JNT, Nj^Hs)' The large 
range in kinematic distances results in a large range in 
most physical properties (see Table [S] for a full statisti- 
cal description of the distribtuion of each physical quan- 
tity). Different scale structures are detected at different 
distances due to the Malmquist bias and the spatial fil- 



tering in the BGPS data, which is a complication that 
will affect all continuum based studies. Based on the 
characteristic mass, radius, and density of cores, clumps, 
and clouds presented by Bergin & Tafalla (2007), we find 
that the BGPS predominately detects cores at distances 
of 1 to 2 kpc, clumps at distances of 2-7 kpc, and clouds 
at distances greater than 7 kpc. Therefore, the effects of 
distance must be carefully accounted for. 

We have studied the Galactic trends in each physical 
property as a function of Galactocentric radius, i?Gai- 
We see a strong peak in the number of BGPS sources 
within the 5 kpc molecular ring, which corresponds to 
the end of the Galactic bar and beginning of the Scutum 
spiral arm. We also see a peak in sources at /Jcai^S 
kpc, which corresponds to the Sagittarius Arm. We have 
found the following trends in physical properties with 

• There are no obvious trends among the inner 
Galaxy sources. Sources in Gem OBI have a higher 
(Tk) than sources within the inner Galaxy, al- 
though this could be caused by a bias of the Gem 
OBI sources toward later evolutionary stages. 

• NH3 column density decreases with i?Gai with a 
factor of approximately 8 difference between the 
sources closest to the Galactic center and those in 
Gem OBI. 

• Similarly, NH3 abundance also decreases smoothly 
with i?Gai with a factor of 7 difference between the 
inner and outer most BGPS sources. The NH3 
abundance changes with i?Gai at the same rate as 
the nitrogen abundance, suggesting the decrease in 
NH3 abundance is caused by less nitrogen to form 
NH3 toward the outer Galaxy, although changes in 
dust properties (e.g. gas-to-dust ratio, grain size, 
grain temperature) likely also affect the measured 
NH3 abundance. 

• When comparing sources at similar kinematic dis- 
tances, there is little variation in physical proper- 
ties across the Galaxy regardless of enviornment 
(inside versus outside the molecular ring, for ex- 
ample) . This is in contrast to the molecular clouds 
seen in the GRS, which were found to be larger, 
more massive, warmer, and have higher column 
densities within the molecular ring (Rathborne et 
al. 2009). This suggests that environment may af- 
fect the large scale molecular clouds while interior 
regions of dense gas exhibit similar properties re- 
gardless of environment. 

We have characterized the star formation activity of 
the BGPS sources based on the presence of an associated 
mid-IR source as described in Dunham et al. (2011). We 
find that 67% of the BGPS sources in this study have an 
associated mid-IR source. BGPS sources associated with 
either an Extended Green Object (EGO; Cyganowski et 
al. 2008) or a Red MSX Source (RMS; Hoare et al. 2004; 
Urquhart et al. 2008) have a higher mean Tk, cryisr, and 
-^NH3 than sources without any mid-IR indication of 
star formation activity. 

We also compared the BGPS sources to criteria sug- 
gested for both efficient star formation and for massive 
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star formation. We find that 46.3% have average surface 
densities above the threshold for efficient star formation 
(Heiderman et al. 2010; Lada et al. 2010). If we con- 
sider the peak column density, we find 86.8% of the full 
sample have some part with surface densities above the 
threshold for efficient star formation. Similarly, we have 
compared the BGPS sources to the criterion for massive 
star formation presented by KaufFmann et al. (2010), and 
found that 47.6% of the full sample have masses above 
their criterion. Thus, roughly half of the full sample sat- 
isfies both criteria and should be forming stars with high 
efficiency, while smaller parts of 86.8% of the full sam- 
ple may be forming stars efficiently. This is consistent 
with the 67% of the full sample that has correspond- 
ing mid-IR sources. The BGPS sources that satisfy the 
Kauffnian et al. (2010) criterion for massive star forma- 
tion are, on average, larger, more massive, and less dense 
than seen in previous large-scale studies of massive star- 
forming regions that were based on tracers of massive 
star formation (i.e. Plume et al. 1992, 1997; Shirley et 
al. 2003; Wu et al. 2010). 
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TABLE 4 
Observed NH3 Properties 



ID RA 


Dec 


Vlsr 


'^Vlsr T^b(l,l) VF(1,1) 


rn,b(2,2)b 


W^(2,2)b.c 


T„ib(3,3)b.<^ 


W'(3,3)'''<= 


rn,b(4,4)b.<= 


W^(4,4)b,c 


Number" (J2000) 


(J2000) 


(km s-i) 


(kms-i) (K) (Kkms-l) 


(K) 


(K km s-l) 


(K) 


(K km s-i) 


(K) 


(K km s-l) 



1307 


18 


02 


30.0 


-22 


28 06.9 


1309 


18 


01 


00.5 


-22 


15 48.5 


1310 


18 


02 


47.8 


-22 


26 03.3 


1311 


18 


03 


15.0 


-22 


27 09.7 


1313 


18 


02 


45.3 


-22 


21 46.8 


1314 


18 


03 


11.9 


-22 


24 32.5 


1316 


18 


03 


31.5 


-22 


26 48.4 


1317 


18 


03 


26.8 


-22 


25 16.1 


1320 


18 


02 


46.6 


-22 


12 17.7 


1322 


18 


03 


23.7 


-22 


07 45.9 



32.63 (0.03) 
8.16 (0.06) 
122.87 (0.03) 
153.30 (0.01) 
-0.71 (0.04) 
154.02 (0.01) 
152.52 (0.00) 
154.55 (0.02) 
18.11 (0.03) 
0.26 (0.00) 



0.94 (0.03) 
0.94 (0.06) 
0.60 (0.02) 
0.80 (0.01) 
0.77 (0.04) 
0.97 (0.01) 
0.63 (0.00) 
1.01 (0.02) 
0.76 (0.03) 
0.42 (0.00) 



0.85 (0.10) 
0.39 (0.09) 
0.75 (0.09) 
2.54 (0.10) 
0.63 (0.09) 
3.38 (0.11) 
3.93 (0.10) 
1.26 (0.10) 
0.75 (0.09) 
3.00 (0.09) 



6.08 (0.20) 
3.21 (0.18) 
4.08 (0.18) 
19.31 (0.21) 
3.94 (0.19) 
28.53 (0.23) 
24.58 (0.20) 
10.45 (0.21) 
4.24 (0.18) 
15.52 (0.16) 



< 0.35 

< 0.35 

< 0.35 
1.38 (0.09) 

< 0.35 

2.07 (0.09) 
2.09 (0.09) 
0.48 (0.09) 

< 0.34 

1.08 (0.09) 



< 0.97 

< 0.97 

< 0.86 
5.88 (0.19) 

< 0.92 
10.59 (0.20) 
7.47 (0.18) 
2.49 (0.20) 

< 0.90 
2.12 (0.15) 



< 0.39 

< 0.39 

< 0.38 
1.02 (0.10) 

< 0.37 
1.60 (0.10) 
1.19 (0.09) 

< 0.38 

< 0.38 

< 0.38 



< 0.95 

< 0.98 

< 0.83 
7.09 (0.19) 

< 0.88 
12.32 (0.20) 
6.53 (0.17) 

< 0.96 

< 0.89 

< 0.74 



Note. — Errors are given in parentheses. The full table is available in the online journal. 
" B, C, D, and E denote multiple ammonia pointings that fall within a single 1.1 mm source. 

^ Upper-limits are T^(, < 4(t and W < 5(tAi;\/^, where cr is the rms noise, Av is the width of a single channel in velocity, and N is the number of pixels over which the average RMS was 
calculated. 

^ • • • denotes sources that were not observed in the given transition. 
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TABLE 5 
Derived Gas Properties 



ID 






7^ 
^ ex 


(X 


a NT 


X12 W 


A/T 1 1 1 + i T-il o 
iViUibipit; 


Number^ 


r(l,l) 


(K) 


(K) 


(km s-1) 


(km 


>s-i) 


maser?*^ 


Comp?d 


1307 


4.31 


< 11.90 (0.47) 


3.29 (0.13) 


0.20(0.002) 


0.93 


(0.06) 


Y 




1309 


2.70 


< 28.59 (1.17) 


3.01 (0.13) 


0.32(0.004) 


0.93 


(0.11) 


Y 




1310 


4.31 


< 12.61 (0.46) 


3.27 (0.11) 


0.21(0.002) 


0.59 


(0.05) 


N 




1311 


6.06 


24.29 (0.19) 


4.03 (0.04) 


0.29(0.001) 


0.80 


(0.02) 


Y 


Vlsr 


1313 


2.79 


< 14.00 (0.75) 


3.24 (0.17) 


0.22(0.003) 


0.77 


(0.08) 


N 




1314 


5.12 


27.44 (0.13) 


4.60 (0.03) 


0.31(0.000) 


0.96 


(0.01) 


Y 


Tk 


1316 


5.15 


18.01 (0.10) 


5.51 (0.04) 


0.25(0.000) 


0.63 


(0.01) 


Y 




1317 


3.62 


14.40 (0.23) 


3.72 (0.06) 


0.22(0.001) 


1.01 


(0.04) 


N 




1320 


2.14 


< 15.76 (0.68) 


3.42 (0.15) 


0.23(0.002) 


0.76 


(0.07) 


Y 




1322 


8.34 


12.38 (0.12) 


4.79 (0.05) 


0.21(0.000) 


0.41 


(0.01) 


N 


Tk 



Note. — Errors arc given in parentheses. The full table is available in the online journal. 

^ B, C, D, and E denote multiple ammonia pointings that fall within a single 1.1 mm source. 
^ Upper-limits to Tj^j^j^ are due to a non-detection of the NH3(2,2) transition. 
^ Y denotes detected, N denotes not detected. 

*^ Tk denotes additional warm temperature component might provide a better fit, and visr denotes multiple 
velocity components required. 



TABLE 6 
Kinematic Distances 



ID IRDC 


HISA KDA 


kinematic 




Number flag^ 


flagi flag2 


distance (kpc) 


(kpc) 



1307 


n 


n 


f 


i2.4(l'^;g) 


4.1 


1309 


n 


m 


f 




6.6 


1310 


m 


m 


n 


6.9(i;!;is;) 


1.7 


1311 


n 


m 


f 




1.4 


1313 


y 


y 


n 




8.5 


1314 


y 


n 


n 


7.3(i;!;?,?,) 


1.4 


1316 


y 


m 


n 


7 O(- + ().0!K 


1.4 


1317 


y 


m 


n 


7.3(iS;;!!«) 


1.4 


1320 


n 


n 


f 


i3.7(i;;;«^) 


5.3 


1322 


y 


y 




0.1(_}-f) 


8.2 



Note. — Errors are given in parentheses. The full table is 
available in the online journal. 

^ y denotes positive association, n denotes no association, and m 
denotes questionable association. 

^ n denotes near kinematic distance, f denotes far kinematic dis- 
tance, and t denotes tangent distance. 



TABLE 7 
Derived Masses and Densities 



ID 

Number 


Radius 
(pc) 


Miso(120") 


Miso(int) 


Mvir(int) 


np(int) 
(lO^cm-^) 


(lO^cm-^) 


E(int) 
(10-2 E cm-2) 


^NH3 
(lOl* cm-2) 


K7 

(10^1 cm-2) 


(1022"^m-2) 


(10-8) 


1307 


2.1 


o o "1 / r» r\ \ 

/o2i(Dyoj 


2340(580) 


2169(100) 


0.99(0.33) 




3.44(1.01 ) 


10.07(0.93) 


7.35(2.15) 


1.19(0.18) 


13.69(4.20) 


1309 


2.1 


1526(580) 


1632(590) 


2131(210) 


0.74(0.45) 


0.37(0.07) 


2.50(1.22) 


8.04(2.56) 


5.35(2.60) 


0.62(0.07) 


15.01(8.72) 


1310 


1.6 


1163(240) 


985(230) 


665(28) 


0.97(0.24) 


0.76(0.15) 


2.54(0.62) 


6.47(0.76) 


5.43(1.32) 


1.00(0.23) 


11.93(3.23) 


1311 


1.2 


952(150) 


904(120) 


879(11) 


2.33(0.34) 


1.33(0.27) 


4.42(0.64) 


18.41(0.50) 


9.46(1.33) 


0.81(0.09) 


19.45(2.78) 


1313 


0.0 


0(100) 


0(34) 


10(490) 


230.61(81966.01) 


0.80(0.16) 


5.46(0.74) 


5.56(1.07) 


11.67(3278.20) 


0.90(0.21) 


4.76(1336.30) 


1314 


1.2 


932(100) 


1013(90) 


1280(17) 


2.54(0.29) 


2.11(0.42) 


4.87(0.51) 


22.15(0.40) 


10.42(1.06) 


1.68(0.09) 


21.26(2.20) 


1316 


0.7 


745(130) 


719(87) 


302(4) 


10.77(1.53) 


3.79(0.76) 


11.39(1.61) 


14.40(0.22) 


24.37(3.19) 


1.54(0.12) 


5.91(0.78) 


1317 


1.2 


1238(190) 


947(150) 


1474(32) 


2.04(0.37) 


1.46(0.29) 


4.12(0.71) 


10.62(0.61) 


8.81(1.50) 


1.18(0.17) 


12.05(2.17) 


1320 


3.0 


3019(980) 


2923(930) 


2002(140) 


0.45(0.21) 


1.13(0.23) 


2.20(0.86) 


4.53(0.91) 


4.71(1.84) 


0.82(0.18) 


9.61(4.22) 


1322 


0.0 


0(63) 


0(56) 


1(240) 


799.73(929799.74) 


2.44(0.49) 


9.91(9102.62) 


12.08(0.22) 


21.19(19466.69) 


1.77(0.24) 


5.70(5235.07) 



Note. — Errors are given in parentheses. The full table is available in the online journal. 
^ li^x cannot be calculated for sources where we have set TK=Tex. 



TABLE 8 
Statistical Summary 



Pull Sample Tk subsample 



Property 


units 


mimmum 


mean" 


median 


maximum 


mmimum 


mean" 


median 


maximum 


n 


KpC 


u. ± 


7 /I /''^ 


D.O 


10. u 


U. 1 


D.Dl^o. 1 ) 


A 7 




^Gal 


kpc 


1 4 


O . LyL. ^ j 




O.O 


1 4 




A 7 








0.01 


1.20(0.77) 


1.04 


4.79 


0.01 


1.28(0.77) 


0.93 


5.65 




pc 


0.00 


0.79(0.49) 


0.68 


2.75 


0.01 


0.84(0.49) 


0.67 


2.88 




pc 


0.01 


1.87(1.36) 


1.51 


8.41 


0.01 


2.02(1.36) 


1.40 


8.55 


log(S,.(120")) 


Jy 


-0.55 


0.06(0.32) 


0.00 


1.35 


-0.50 


0.22(0.34) 


0.15 


1.35 


log(S„(lnt)) 


Jy 


-0.88 


0.06(0.40) 


0.01 


1.48 


-0.54 


0.24(0.41) 


0.16 


1.48 




K 


0.25 


1.27(0.98) 


0.92 


7.14 


0.38 


1.93(1.12) 


1.70 


7.14 


Tn,b{2,2) 


K 


0.17 


0.59(0.54) 


0.41 


4.73 


0.24 


0.91(0.68) 


0.73 


4.73 


Tex(l,l) 


K 


2.82 


4.07(2.15) 


3.58 


24.81 


2.99 


4.45(2.11) 


4.12 


24.12 




K 


5.00 


15.57(5.02) 


14.24 


59.41 


11.95 


17.39(5.49) 


15.70 


59.41 


(TK-Tex) 


K 


0.00 


11.50(5.05) 


10.38 


56.17 


0.00 


12.94(5.59) 


11.55 


56.17 


Vlsr 


km 


-0.71 


63.7( 32.9) 


60.6 


154.55 


0.26 


59.6( 33.3) 


49.0 


154.5 




km 


0.10 


0.76(0.49) 


0.64 


4.42 


0.24 


0.81(0.53) 


0.69 


4.42 


a 


km 


0.13 


0.23(0.03) 


0.22 


0.46 


0.20 


0.24(0.03) 


0.23 


0.46 


<JJ\TT 


km 


0.08 


0.75(0.49) 


0.63 


4.41 


0.23 


0.81(0.53) 


0.68 


4.41 


Mach {(Jvisr 1 ^) 




0.59 


3.23(1.87) 


2.83 


19.38 


1.12 


3.24(1.73) 


2.85 


15.47 






0.01 


3.44(1.73) 


3.37 


12.22 


0.02 


4.07(1.79) 


3.87 


12.22 


log(Miso(120")) 


M© 


-1.02 


2.95(0.61) 


3.02 


4.40 


-1.02 


2.93(0.62) 


2.94 


4.00 


log(Afiso(int)) 


Mq 


-1.07 


2.94(0.66) 


2.98 


4.39 


-1.07 


2.95(0.67) 


2.98 


4.31 


log(Mvir(int)) 


M0 


-0.07 


2.85(0.66) 


2.88 


4.72 


0.18 


2.87(0.64) 


2.90 


4.69 


log(Mvir/Miso(int)) 




-1.38 


-0.09(0.45) 


-0.13 


1.97 


-0.91 


-0.08(0.42) 


-0.13 


1.97 


log(n(int)) 


cm~^ 


2.07 


3.12( 0.56) 


3.03 


5.90 


2.15 


3.27( 0.58) 


3.22 


5.90 


S(mt) 


g cm^"^ 


0.006 


0.037(0.030) 


0.029 


0.236 


0.010 


0.046(0.033) 


0.036 


0.231 


log(rtei:) 


cm~^ 


2.08 


3.10( 0.35) 


3.09 


4.88 


2.58 


3.28( 0.26) 


3.29 


4.03 


log(rtei:/n(iiit)) 




-2.51 


-0.02(0.58) 


0.06 


2.15 


-2.51 


0.01(0.55) 


0.09 


1.17 


log(Af]vH3 ) 


cm~^ 


13.24 


14.75( 0.37) 


14.76 


15.96 


13.73 


14.95( 0.29) 


14.92 


15.96 


log(^r/) 


cm~^ 


21.08 


21.81( 0.27) 


21.79 


22.70 


21.32 


21.91( 0.27) 


21.89 


22.69 


logCiVlf^"™) 


cm~-^ 


21.30 


22.09( 0.27) 


22.06 


23.29 


21.55 


22.22( 0.28) 


22.18 


23.29 


log(W;f3/iV|p^'"") 




-8.63 


-7.34( 0.30) 


-7.33 


-6.34 


-8.36 


-7.27( 0.26) 


-7.27 


-6.56 



^ Standard deviation is given in parentheses. 
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TABLE 9 

Trends with Star Formation Activity 



Property Probability standard 
(units) Group minimum mean deviation median maximum 

(km s-l) 



(K) 



NH3 



(IQi 






0.1 


0.7 


0.4 


0.6 


3.6 


1 


0.2 


0.7 


0.5 


0.6 


4.4 


2 


0.2 


0.7 


0.5 


0.7 


3.8 


3 


0.3 


1.1 


0.5 


1.0 


3.1 





6.6 


13.9 


3.0 


13.5 


26.6 


1 


5.0 


14.6 


3.8 


13.8 


33.7 


2 


10.3 


15.8 


3.9 


14.6 


28.6 


3 


11.1 


22.7 


8.2 


21.1 


59.4 





0.2 


5.6 


3.8 


4.6 


21.1 


1 


0.2 


7.1 


5.5 


5.8 


40.2 


2 


0.2 


8.0 


6.6 


6.1 


40.0 


3 


2.0 


16.3 


16.9 


10.4 


91.2 



